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SUMMARY 

A  summary  review  is  made  of  the  theory  of  aerodynamic  sound  in  low  Mach  number  flows, 
with  particular  emphasis  on  procedures  for  estimating  the  influence  of  solid  boundaries  on 
sound  production.  Special  consideration  is  given  to  the  production  of  sound  by  nominally 
steady  flow  over  a  lifting  or  control  surface.  Four  categories  of  interactions  are  discussed 
involving  (i)  a  large,  nominally  flat  surface,  (ii)  a  large  surface  with  a  trailing  edge,  (iii)  a 
lifting  surface  of  compact  chord  and  finite  span,  and  (iv)  a  trailing  edge  with  a  streamwise 
slot  or  part-span  flap. 

The  relative  inefficiency  with  which  hydrodynamic  energy  is  converted  into  sound  means 
that  it  is  not  usually  possible  to  make  confident  noise  predictions  by  direct  numerical  simula¬ 
tion  of  a  flow.  However,  numerical  codes  can  accurately  predict  mean  properties  of  complex 
fluid  structure  interactions,  and  the  results  may  be  incorporated  into  analytical  models  de¬ 
scribed  in  this  report  to  express  the  sound  and  vibration  generated  by  an  interaction  in  terms 
of  those  mean  properties.  Mean  flow  numerical  computations  also  permit  ‘noisy’  large  scale 
features  (such  as  trailing  vortices)  to  be  identified,  whose  potential  for  noise  generation  can 
be  investigated  analytically  by  the  methods  described  herein. 
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1.  INTRODUCTION 

The  sound  generated  by  unsteady  flow  is  called  aerodynamic  sound  [1,  2].  Most  flows  of 
technological  interest  are  of  high  Reynolds  number  and  turbulent,  and  the  acoustic  radiation 
is  a  very  small  by-product  of  the  motion.  The  turbulence  is  usually  produced  by  fluid  motion 
relative  to  solid  boundaries  or  by  the  instability  of  free  shear  layers  separating  a  high  speed 
flow  (such  as  a  jet)  from  a  stationary  environment.  The  theory  of  aerodynamic  sound  is 
usually  concerned  with  the  mechanism  of  sound  production  by  a  prescribed  turbulent  flow 
as  opposed  to  the  processes  by  which  the  flow  is  generated.  In  many  instances,  however,  the 
backreaction  of  the  sound  on  the  flow  can  be  important,  and  must  be  properly  accounted 
for  in  estimating  the  radiation. 

The  aerodynamic  sound  problem  in  “free  space”  reduces  to  the  study  of  mechanisms  that 
convert  kinetic  energy  of  free  rotational  (“vortical”)  motions  into  acoustic  waves  involving 
longitudinal  vibrations  of  fluid  particles.  There  are  two  principal  source  types:  a  quadrupole, 
whose  strength  is  determined  by  the  unsteady  Reynolds  stress,  and  a  dipole,  which  is  im¬ 
portant  when  mean  mass  density  variations  occur  within  the  source  region. 

The  quadrupole  sound  generated  by  turbulent  flow  near  a  fixed,  rigid  surface,  is  augmented 
by  radiation  from  a  distribution  of  surface  dipoles  whose  strength  is  the  force  per  unit  sur¬ 
face  area  exerted  on  the  fluid.  If  the  surface  is  in  accelerated  motion  there  are  additional 
dipoles  and  quadrupoles;  neighboring  surfaces  in  relative  motion  experience  “potential  flow” 
interactions  that  also  generate  sound.  At  the  low  Mach  numbers  M  encountered  in  underwa¬ 
ter  applications  (typically  <  0.01),  the  acoustic  efficiency  of  the  surface  sources  exceeds  the 
quadrupole  efficiency  by  a  large  factor  ~  0(1/M2).  The  presence  of  solid  surfaces  within  low 
Mach  number  turbulence  can  therefore  lead  to  substantial  increases  in  aerodynamic  sound 
levels. 

Turbulence  stresses  on  an  elastic  surface  generate  sound,  and  can  also  excite  structural 
vibrations  which  store  a  substantial  amount  of  flow  energy.  The  vibrations  are  ultimately 
dissipated  by  frictional  losses,  but  they  can  contribute  significantly  to  the  radiated  noise, 
because  elastic  waves  are  “scattered”  at  structural  discontinuities,  and  some  of  their  energy 
is  transformed  into  sound.  Flow  generated  sound  accordingly  reaches  the  far  field  via  two 
paths:  directly  from  the  turbulence  sources  and  indirectly  from  possibly  remote  locations 
where  the  scattering  occurs.  In  consequence,  the  effective  acoustic  efficiency  of  the  flow  is 


Section  1 


2 


Introduction 


Report  No.  AM  98-001 


Boston  University,  College  of  Engineering 


often  very  much  larger  than  for  a  geometrically  similar  rigid  surface,  even  when  only  a  small 
fraction  of  the  structural  energy  is  scattered  into  sound. 

The  first  part  of  this  report  (§2)  consists  of  a  review  of  the  theory  of  aerodynamic 
sound  for  low  Mach  number  flows,  including  methods  for  estimating  the  contributions  from 
solid  boundaries.  Subsequent  sections  deal  with  the  prediction  of  sound  generated  by  high 
Reynolds  number  flow  over  typical  structural  elements  of  a  complex,  underwater  control 
surface.  The  flow  induced  sound  and  vibration  of  a  submerged  vehicle  is  attributable  to  the 
instability  of  boundary  layers,  the  shedding  of  larger  scale  vortex  structures  from  appendages, 
and  to  localized  sources  formed  when  these  flow  inhomogeneities  interact  with  or  impinge  on 
geometrical  and  structural  discontinuities,  such  as  joints,  corners,  fixed  and  articulated  lift¬ 
ing  surfaces,  shrouds,  etc.  With  few  exceptions,  and  because  of  the  large  number  of  degrees 
of  freedom,  and  the  relative  inefficiency  with  which  hydrodynamic  energy  is  converted  into 
sound,  it  is  not  possible  to  make  confident  noise  predictions  by  direct  numerical  calculation 
of  the  flow. 

However,  because  of  the  increasing  ability  of  numerical  codes  to  simulate  accurately  the 
principal  mean  properties  of  fluid  structure  interactions,  the  results  of  such  calculations  can 
be  incorporated  into  analytical  prediction  schemes  to  estimate  the  sound  and  vibration  gen¬ 
erated  by  fluid-structure  interactions  in  terms  of  computed  mean  hydrodynamic  properties. 
For  example,  numerical  solutions  of  the  Reynolds-Averaged  Navier-Stokes  (RANS)  equa¬ 
tions,  which  make  use  of  a  simplified  turbulence  model,  enable  the  mean  flow  field  about 
geometrically  complex  structures  to  be  evaluated  in  detail.  Such  solutions  can  be  used  to 
calculate  characteristic  ‘wavelengths’,  wavenumbers  and  frequencies  for  insertion  into  empir¬ 
ical  models  of  the  unsteady  surface  pressure  and  shear  stress  spectra  that  are  then  used  to 
estimate  sound  and  vibration  levels  produced  by  boundary  layer  forcing  of  an  elastic  sub¬ 
structure.  These  methods  also  permit  ‘noisy’  large  scale  features  (such  as  trailing  vortices)  to 
be  identified,  whose  unsteady  interaction  with  a  structural  appendage  or  surface  irregularity 
can  be  investigated  by  appeal  to  analytical  models  involving  the  perturbation  of  the  large 
scale  structure  by  a  turbulence  field,  or  from  some  known  instability  characteristic  of  the 
motion. 

The  purpose  of  this  report  is  to  review  and  extend  where  necessary  the  available  analytical 
models  for  low  Mach  number  flows  interacting  with  lifting  surfaces.  This  is  done  by  consid¬ 
ering  four  different  classes  of  flow  interactions  with:  (i)  large,  nominally  flat  surfaces  (§3), 
(ii)  a  large  surface  with  a  trailing  edge  (§4),  (iii)  lifting  surfaces  of  finite  chord  and  span  (§5), 
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and  (iv)  trailing  edges  with  streamwise  slots  and  part-span  flaps  (§6).  Well  known  methods 
formulae  are  reviewed  without  extensive  discussion,  but  results  not  available  elsewhere 
in  the  literature  are  treated  in  more  detail  (especially  in  §§4-6). 
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2.  THEORY  OF  AERODYNAMIC  SOUND 
§2.1  Lighthill’s  theory 

The  acoustic  analogy.  Lighthill  [1]  reformulated  the  Navier-Stokes  equation  into  an  exact, 
inhomogeneous  wave  equation  whose  source  terms  are  important  only  within  the  turbulent 
(vortical)  region.  Sound  is  expected  to  be  such  a  very  small  component  of  the  flow  that,  once 
generated,  its  back-reaction  on  the  source  flow  is  usually  negligible.  In  a  first  approximation 
the  motion  in  the  source  region  may  therefore  be  determined  by  neglecting  the  production 
and  propagation  of  the  sound.  This  would  be  inappropriate  if  the  Mach  number  M  is  large 
enough  for  compressibility  to  be  important  in  the  source  flow,  when  the  source  flow  is  coupled 
to  a  resonator,  such  as  an  organ  pipe,  or  when  bubbles  are  present  in  the  case  of  liquids. 
However,  there  are  many  technologically  important  flows  where  M  is  sufficiently  small  that 
the  hypothesis  is  obviously  correct,  and  where  the  theory  leads  to  unambiguous  predictions 
of  the  sound. 

Consider  the  sound  generated  by  a  finite  region  of  rotational  flow  in  an  unbounded  fluid 
at  rest  at  infinity.  Let  us  compare  the  equations  for  the  density  fluctuations  in  the  real 
flow  with  those  for  an  ideal,  linear  acoustic  medium  which  coincides  with  the  real  fluid  at 
large  distances  from  the  sources.  The  difference  between  these  equations  will  be  shown  to  be 
equivalent  to  a  distribution  of  sources  in  the  ideal  acoustic  medium,  whose  radiation  field  is 
the  same  as  that  in  the  real  flow  and  may  therefore  be  calculated  by  the  methods  of  linear 
acoustics.  To  do  this,  body  forces  are  neglected,  and  the  momentum  equation  is  written  in 
the  form 

d(pvi)/dt  =  -ditij/dxj,  (2.1.1) 

where  7Ty  is  the  momentum  flux  tensor 

7Ty  =  pViVj  +  (p  -  p0)Sij  -  Oij,  (2.1.2) 

p0  is  the  uniform  pressure  at  infinity,  and  ay  is  the  viscous  stress  tensor  (other  standard 
symbols  are  defined  on  page  92). 

By  integrating  over  a  fixed  region  V,  it  can  be  seen  that  (2.1.1)  equates  the  rate  of  change 
of  momentum  in  V  to  the  action  of  the  pressure  and  viscous  stresses  on  its  boundary  and 
to  the  convection  of  momentum  across  the  boundary  at  a  rate  determined  by  the  Reynolds 
stress  pv^j.  In  an  ideal,  linear  acoustic  medium,  momentum  transfer  is  produced  solely 
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by  the  pressure,  i.e.,  tt^  -*  tt°  =  (p  -  p0)<J#  =  c*(p  -  p0)6tj,  where  p0  and  c0  are  the  mean 
density  and  sound  speed.  The  equation  of  linear  acoustics  for  the  perturbation  density  p-p0 
is  therefore 

(d2/clat2  -  v2)  [«*(/>  -  ft,)]  =  o,  (2.1.3) 

and  is  obtained  by  eliminating  pv(  between  (2.1.1)  (with  7r replaced  by  7r?)  and  the  conti¬ 
nuity  equation 

dp/dt  +  div(pv)  =  0.  (2.1.4) 

In  the  absence  of  externally  applied  forces  or  moving  boundaries,  equation  (2.1.3)  has  only 
the  trivial  solution  p  -  p0  =  0,  since  the  radiation  condition  ensures  that  sound  waves  cannot 
enter  from  infinity. 

The  sound  generated  in  the  real  fluid  may  now  be  seen  to  be  exactly  equivalent  to  that 
produced  in  an  ideal,  stationary  acoustic  medium  that  is  forced  by  the  stress  distribution 
Tij  =  7 Tij  -  7 r°j.  Tij  is  the  Lighthill  stress  tensor 

Tij  =  pviVj  +  ((p  -  Po)  -  C20{p  -  po ))  Sij  -  Gij.  (2.1.5) 

The  Reynolds  stress  pViVj  is  nonlinear,  and  is  significant  only  within  the  rotational  source 
region.  The  second  term  is  the  excess  of  momentum  transfer  by  the  pressure  over  that  in 
an  ideal  (“linear”)  fluid  of  density  p0  and  sound  speed  c0.  This  is  caused  by  wave  amplitude 
nonlinearity,  and  by  mean  density  variations  in  the  source  flow.  The  viscous  stress  tensor 
Gij  is  linear  in  the  perturbation  quantities,  and  properly  accounts  for  the  attenuation  of  the 
sound;  in  most  applications  the  Reynolds  number  in  the  source  region  is  sufficiently  large 
that  g^  can  be  neglected,  and  attenuation  in  the  radiation  zone  is  usually  ignored  in  a  first 
approximation. 

Lighthill’s  acoustic  analogy  equation  for  the  production  of  aerodynamic  sound  is  ob¬ 
tained  by  first  re-writing  (2.1.1)  as  the  momentum  equation  for  an  ideal,  stationary  fluid  of 
density  p0  and  sound  speed  c0  subject  to  the  externally  applied  stress  T^: 


d(pvi)/dt  +  d  (c20{p  -  po ))  /dxi  =  -dTij/dxj. 


Elimination  of  the  momentum  density  pV{  between  this  and  the  continuity  equation  (2.1.4) 
yields  Lighthill’s  equation,  which  is  the  exact,  nonlinear  counterpart  of  (2.1.3): 


\c2dt 2 


dxidxj 


(2.1.6) 
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The  fluid  mechanical  problem  of  calculating  the  aerodynamic  sound  is  therefore  formally 
equivalent  to  solving  this  equation  for  the  radiation  into  a  stationary,  ideal  fluid  produced 
by  a  distribution  of  quadrupole  sources  whose  strength  per  unit  volume  is  the  Lighthill  stress 
tensor  7y. 

The  solution  with  outgoing  wave  behavior  in  an  unbounded  fluid  is  given  by  the  retarded 
potential  integral 


clip-  Po ) 


d2  [Tjj(y,t-  |x-y|/cp) 

dxidxj  J  47t|x  —  y| 


(2.1.7) 


where  clip  ~  Po)  P~Po  as  |x|  -»  oo  in  the  linearly  disturbed  fluid  outside  the  source  flow. 
This  solution  should  strictly  be  regarded  as  an  alternative,  integral  equation  formulation 
of  the  Navier-Stokes  equation,  that  happens  to  provides  a  useful  representation  of  the  far 
field  sound  when  Ty  is  known.  Ty  accounts  not  only  for  the  generation  of  sound,  but  also 
for  self-modulation  due  to  acoustic  nonlinearity,  convection  by  the  flow,  refraction  due  to 
sound  speed  variations,  and  attenuation  due  to  thermal  and  viscous  actions.  Nonlinear 
effects  on  propagation  and  dissipation  are  usually  sufficiently  weak  to  be  neglected  within 
the  source  region,  although  they  may  affect  propagation  to  a  distant  observer.  Convection 
and  refraction  of  sound  within  and  near  the  source  flow  can  be  important,  for  example  when 
the  sources  are  contained  in  a  turbulent  shear  layer,  or  are  adjacent  to  a  large,  quiescent 
region  of  fluid  whose  mean  thermodynamic  properties  differ  from  those  in  the  radiation 
zone.  Effects  of  this  kind  are  accounted  for  by  contributions  to  7y.  that  are  linear  in  the 
perturbation  quantities  relative  to  a  mean  background  flow  [2]. 

Thus,  with  the  exception  of  flows  amenable  to  special  treatment,  the  practical  utility  of 
LighthilPs  equation  rests  on  the  hypothesis  that  all  of  these  effects,  which  actually  depend 
on  the  compressibility  of  the  source  flow,  can  be  ignored,  and  that  adequate  predictions  of 
the  aerodynamic  sound  are  obtained  by  taking  for  7y  an  estimate  based  on  the  equations 
of  motion  of  an  incompressible  fluid.  This  approximation  is  likely  to  be  acceptable  when 
M2  <C  1,  and  when  the  wavelength  of  the  sound  is  much  larger  than  the  dimension  of  the 
source  region.  The  remainder  of  this  section  is  devoted  to  a  consideration  of  such  cases. 

Aerodynamic  sound  generated  by  low  Mach  number  turbulence  of  uniform  mean 
density.  When  the  mean  density  and  sound  speed  are  uniform,  the  variations  in  p  produced 
by  low  Mach  number,  high  Reynolds  number  velocity  fluctuations  are  of  order  p0M 2,  and 
pViVj  «  p0ViVj  with  a  relative  error  ~  0(M2)  <C  1.  Similarly,  p  -  pQ  -  clip  —  p0)  ~ 
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(p-p0)(  1  -  c20/c2)  ~  0(p0v2M2).  Thus,  Tij  «  p0ViVj ,  when  viscous  stresses  are  neglected, 
and  the  solution  (2.1.7)  becomes 


P(x,t) 


d2  r  poVjVj(y2 1  -  jx  -  yj/cp)  rf3 

foi&Cj  J  47r|x  -  y| 


«  /  /W;(y>  t  -  |x  -  y|/co)  d3y,  |x|  ->  00,  (2.1.8) 

where  p(x,t)  =  c20(p  -  p0)  is  the  perturbation  pressure  in  the  far  field.  Quantitative  pre¬ 
dictions  can  be  made  from  this  equation  provided  the  behavior  of  the  Reynolds  stress  is 
known. 

The  order  of  magnitude  of  p  can  be  estimated  in  terms  of  the  characteristic  velocity 
v  and  length  scale  £  (of  the  energy  containing  eddies)  in  the  source  region.  £  is  deter¬ 
mined  by  the  scale  of  the  mechanism  responsible  for  turbulence  production,  such  as  the 
width  of  a  jet  mixing  layer.  Fluctuations  in  vivi  occurring  in  different  regions  of  the  tur¬ 
bulent  flow  separated  by  distances  >  0(£)  will  tend  to  be  statistically  independent,  and 
the  sound  may  therefore  be  considered  to  be  generated  by  a  collection  of  V0/£3  inde¬ 
pendent  eddies,  where  V0  is  the  volume  occupied  by  the  turbulence.  The  dominant  fre¬ 
quency  of  the  motion  ~  u/^,  so  that  the  wavelength  of  the  radiated  sound  ~  £/M  £ 

(M  ~  v/c0),  and  each  eddy  is  therefore  acoustically  compact.  Hence,  the  acoustic  pressure 
generated  by  a  single  eddy  is  of  order  p  ~  (£/\x\ )p0v2M2,  and  the  acoustic  power  it  radiates 
~  47t|x| 2p2/p0c0  ~  £2p0v*/c50  =  £2p0vzMb.  This  is  Lighthill’s  “eighth  power”  law.  The  total 
acoustic  power  is  na  «  (' V0/£3){£2p0vzM 5)  =  p0vzMhV0/£.  Dimensional  arguments  and  ex¬ 
periment  [3,  4]  indicate  that  the  rate  of  decay  of  the  turbulence  kinetic  energy  ~  0(V0p0vz/£). 
In  a  statistically  steady  state  this  must  equal  the  rate  IT0,  say,  at  which  energy  is  supplied 
to  the  flow  by  the  action  of  external  forces.  The  efficiency  na/no  with  which  this  energy  is 
converted  into  sound  is  therefore  proportional  to  M5,  confirming  Lighthill’s  hypothesis  that 
the  flow  generated  sound  is  an  infinitesimal  by-product  of  the  motion. 

Aerodynamic  sound  generated  by  low  Mach  number  turbulence  of  variable  mean  density.  When 
the  mean  density  in  the  source  region  is  not  constant,  the  Reynolds  stress  quadrupoles  are  augmented  by 
a  distribution  of  dipoles,  the  dipole  strength  being  proportional  to  the  hydrodynamic  force  experienced  by 
a  fluid  particle  of  density  p  relative  to  the  force  the  same  particle  would  experience  had  its  density  been 
equal  to  pQ.  This  dipole  is  “hidden”  within  the  Lighthill  quadrupole  Ty,  and  is  estimated  by  first  writing 
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the  solution  (2.1.7)  in  the  form 


p(x,t) 


1  d2 
Aitc2  dt2 


J  [pu.Vj] 


(Xj  -  yi)(xj  -  yj)<Py  ,  1  f  \p~Po  _(f)_ 

J  dt2  [ 


|x  -  y|3 


47T 


( P  -  Po ) 


d3y 

|x-y|’ 


|x  -  y|  ->  oo 
(2.1.9) 


where  the  square  brackets  [  ]  indicate  evaluation  at  the  retarded  time  t  -  |x  -  y|/c0. 


The  density  of  a  fluid  particle  changes  by  (p  -  p0)/c 2  when  isentropically  compressed  by  departures  of  p 
from  the  ambient  pressure  p0,  where  c  is  the  local  speed  of  sound.  Let  us  assume  that  the  motion  of  each  fluid 
particle  is  isentropic,  but  that  the  mean  density,  which  will  be  denoted  by  p(x,  t),  varies  from  point  to  point 
in  the  source  region,  where  (p  -  p0)/po  need  not  necessarily  be  small.  This  means  that  (p  -p0)/c2  =  p-p 
and  Dp/Dt  =  0.  The  Reynolds  stress  on  the  right  (2.1.9)  must  be  replaced  by  pviVj,  and  the  integrand  of 
the  second  term  in  (2.1.9)  becomes 


d2  f  P  ~  Po 
dt2  1  c2 


(P  -  Po) 


( P-Po ) 


a2 

dt2 


(p  -  Po)- 


The  first  term  on  the  right  is  already  in  a  form  suitable  for  estimating  the  acoustic  radiation.  However,  the 
second  term  would  yield  zero  when  retarded  time  variations  are  neglected  over  a  coherent  source  region  of 
constant  mass.  To  deal  with  this  the  relation  D(p  —  pa)/Dt  =  0  is  used  to  show  that 


a2 .  ,  _ a_ 

gt2(p  Po)  fa. 


+  dx  dx~  ((P  -  Po)viVj)  +  £  ((p  -  Po)div  v)-£-  (vj(p  -  po) div  v) . 
Zt  X]  3  (2.1.10) 


This  is  further  simplified  by  substitution  from  the  momentum  equation  Dvj/Dt  =  (~l/p)dp/dxj,  and  from 
the  continuity  equation  divv  =  (-1  /pc2)Dp/Dt  «  -D[{p- p0)/pc2]/Dt.  The  amplitude  of  sound  produced 
by  the  final  term  on  the  right  of  (2.1.10)  is  0{M)  relative  to  the  preceding  one  and  is  neglected. 


Inserting  these  results  into  (2.1.9),  and  discarding  the  overbar  on  p,  we  obtain 


p(x,t) 


XiXj  a2  f  r  l 

+-2s_il 

47r|x|  at2  y  L\p0c^ 


d3  y + 


PoXj  d  f  [/_!__  1\  ap‘ 
4nc0\x\2dtJ  [VPo  p)  dyj 


|x|  -4  oo, 


d3  y 


(2.1.11) 


where  quantities  in  square  braces  are  evaluated  at  time  t  —  |x  —  y|/c0.  The  new  terms  in  this  formula,  in 
addition  to  the  Reynolds  stress  radiation  of  (2.1.8),  are  respectively  non-zero  when  the  mean  density  in  the 
source  region  p  ^  p0,  and  when  p0c2  ^  pc2,  and  are  the  sound  fields  of  dipole  and  monopole  sources. 


All  of  these  source  types  occur  in  the  turbulent  mixing  region  of  a  hot  gas  jet  exhausting  into  cold  air:  “hot 
spots”  or  “entropy  inhomogeneities”  behave  as  scattering  centers  at  which  dynamic  pressure  fluctuations  are 
converted  directly  into  sound.  The  dipole  source  strength  is  proportional  to  (1/p  —  l/p0)Vp,  i.e.,  to  the 
difference  between  the  acceleration  of  fluid  of  density  p  in  the  jet  and  that  which  fluid  of  ambient  mean 
density  p0  would  experience  in  the  same  pressure  gradient.  For  an  ideal  gas  with  p  -p0  ~  0(pov2)  in  the 
jet,  the  order  of  magnitude  of  the  dipole  sound  from  an  eddy  of  scale  £  is 


p  ~  (£/|x|)(l  -  p0/p)p0v2M  *  (£/\x\)(AT/T)p0v2 M, 


where  A T/T  ~  (p  -  Po)  I P  is  the  fractional  temperature  difference  between  the  hot  spot  and  its  environment, 
which  can  be  large.  The  sound  power  ~  £2(AT/T)2p0v3M3  exceeds  that  from  a  Reynolds  stress  quadrupole 
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by  a  factor  of  order  (A T/T)2/M2.  Thus,  “entropy  noise”  may  be  an  important  component  of  the  noise  of  a 
hot  gas  jet  at  very  low  Mach  numbers  [5]. 

The  final  term  on  the  right  of  (2.1.11)  is  a  monopole  whose  strength  is  determined  by  the  difference 
between  the  adiabatic  compressibilities  ( K ,  =  1  /pc2)  in  the  source  region  and  in  the  ambient  medium.  In 
an  ideal  gas  Ks  =  1/7 p,  so  that  this  source  is  generally  small,  although  it  may  be  significant  during  the 
turbulent  mixing  of  gases  with  different  values  of  the  specific  heat  ratios  7.  It  can  be  very  important  in 
multiphase  flows  where,  for  example,  the  presence  of  small  air  bubbles  in  water  often  leads  to  an  immense 
increase  in  the  turbulence  generated  noise  [6]. 

§2.2.  Integral  representation  of  aerodynamic  sound 

Most  practical  problems  of  sound  generation  by  flow  involve  moving  boundaries,  moving 
sources  interacting  with  boundaries,  or  turbulence  in  shear  layers  separating  a  quiescent 
medium  from  a  high  speed  flow.  To  apply  Lighthill’s  equation  in  these  circumstances  control 
surfaces  S  are  introduced.  These  may  coincide  with  the  surface  of  a  moving  solid,  or  mark 
a  convenient  interface  between  fluid  regions  of  widely  differing  mean  properties.  A  solution 
is  then  sought  by  imposing  boundary  conditions  on  S,  either  by  first  performing  subsidiary 
calculations  to  determine  the  pressure  or  velocity  on  S  or,  when  S  coincides  with  the  surface 
of  a  solid,  by  application  of  suitable  impedance  conditions. 

Let  /(x,  t)  be  an  indicator  function  that  vanishes  on  S  and  satisfies  /(x,t)  >  0  in  the  fluid 
where  LighthilPs  equation  is  to  be  solved,  and  /(x,t)  <  0  elsewhere,  and  set 

Pij  =  {P~  P0)Sij  ~  Oiy  (2.2.1) 

Multiply  the  momentum  equation  (2.1.1)  by  the  Heaviside  step  function  H(/)  and  rearrange 
into  the  form 

jt  {pvi H(/))  +  ^7  (pViK{f)cl(p  ~  Po))  =  -^7  (H (f)Tij)  +  (pvi{vj  -  vj)  +Pij)^-.(f)  > 

where  v  is  the  velocity  of  S  (so  that  dH(f)/dt  =  —VjdH.(f)/dXj),  and  p0,  c0  are  the  mean 

density  and  sound  speed  in  /  >  0.  To  simplify  the  notation  we  shall  henceforth  write  H 

instead  of  H(/)  when  there  is  no  danger  of  confusion.  The  same  procedure  applied  to  the 

equation  of  continuity,  supplies 

ft  ft  /5H 

—  (H(p  -  po))  +  —  (H pvi)  =  (p(vi  -  Vi)  +  p0Vi)  ^7. 

The  elimination  of  Hpu,  between  these  two  equations  yields  the  differential  Ffowcs  Williams 
-  Hawkings  equation 
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dx(dxj  a*Al  t  j  v’l+p'>'dxi. 


(2.2.2) 


This  equation  is  valid  throughout  the  whole  of  space.  By  using  the  free  space  Green’s 
function 

G(x,  y,  t  -  t)  =  4^“^  “  T  ~  lx  “  y  1/ °o)  (2.2.2) 

to  write  down  the  formal,  outgoing  wave  solution,  it  is  transformed  into  an  integral  equation 
known  as  the 


Ffowcs  Williams  -  Hawkings  equation  [7] 


=  ai/vWN4^-|-/sWH,'-^)+^] 


+ 1  ir,  ^  +  (2-2-4) 

where  quantities  in  the  square  braces  [  ]  are  evaluated  at  the  retarded  time  r  =  t- |x-y|/c0, 
the  surface  integrals  are  over  the  retarded  surface  S(r)  defined  by  /(y,  r)  =  0,  with  the 
surface  element  dS  directed  into  the  region  V(r)  where  /  >  0. 

If  the  control  surface  is  stationary  (v  =  0),  S  is  defined  by  f(y,t)  =  f(y)  —  0,  and  (2.2.4) 
reduces  to 


Curie’s  equation  [8] 
Hc>-a,)  = 


al  [  fc.i  h 

dxidxj  Jv  *■  u-*  47r|x  —  y| 


d_ 

dxi 


dSj(y) 

47t|x  -  y| 


+ 


dSj(y)  } 

4-7r|x  -  y|’ 


(2.2.5) 


which  is  a  generalization  of  Kirchhoff’s  formula  of  linear  acoustics.  When  S  is  stationary 
and  rigid,  this  simplifies  to 


Hc^(p  -  p0) 


f  fT..i  d>y  ±  /  [„;  i  ^(yL. 

dx^Xj  Jv  47t|x  —  y|  dxi  Js  *■  13 *  47r|x  —  y| 


(2.2.6) 


The  surface  integrals  in  these  formulae  may  be  interpreted  as  distributions  of  monopole  or 
dipole  sources.  They  are  associated  with  mass  and  momentum  transfer  across  S,  but  their 
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values  cannot  normally  be  specified  independently.  Indeed,  because  the  acoustic  power  is 
usually  a  mere  fraction  of  the  total  power,  small  errors  made  in  prescribing  the  surface  terms 
can  seriously  impair  the  accuracy  of  predictions  of  the  flow-generated  sound. 


§2.3.  Theory  of  vortex  sound 

The  role  of  vorticity  in  Lighthill’s  theory.  At  low  Mach  numbers  the  velocity  defining 
the  Reynolds  stress  quadrupole  in  Lighthill’s  equation  can  be  determined  by  regarding  the 
source  flow  as  incompressible,  and  using  the  Biot-Savart  law  [9]  to  express  v  in  terms  of  the 
vorticity.  This  suggests  the  existence  of  a  fundamental  relation  between  the  sound  waves, 
which  diverge  from  the  source  region  and  transport  energy  to  the  whole  of  space,  and  the 
vorticity,  which  is  confined  by  Kelvin’s  theorem  [9]  to  a  finite  domain  where  most  of  its 
kinetic  energy  is  ultimately  dissipated  by  the  action  of  viscosity. 


Consider  an  acoustically  compact,  homentropic  vorticity  distribution  u>  of  scale  t  centered 
on  the  coordinate  origin  in  a  medium  of  density  p0  and  sound  speed  c0.  The  velocity  v  = 
u  +  Vp,  where  u  is  defined  in  terms  of  the  vorticity  oj  by  the  Biot-Savart  formula 


/  ,  fv(y,t)d3y 

u(x't)=curl;  toRTTp 


(2.3.1) 


which  is  a  purely  kinematic  relation  between  u)  and  a  velocity  u  that  vanishes  at  infinity. 
Since  /  w(x,t)  dzx  =  0,  it  follows  that  u  ~  0(l/|x|3)  as  |x|  -¥  oo.  Also  div  u  =  0,  so 
that  the  value  of  p  is  determined  by  the  compressibility  of  the  fluid.  In  the  source  region 
P  —  Vo  ~  PoU2,  and  the  characteristic  time  of  the  source  flow  ~  t/u.  Thus,  Dp/Dt  ~  p0u3/£, 
and  the  continuity  equation  S72p  +  Dp/pDt  «  V2<p  +  (1  / p0c20)Dp/ Dt  =  0  implies  that 


Vy?  =  0(uM2),  M  =  u/c0  for  |x|  ~  l. 


(2.3.2) 


We  can  use  the  identity 

d2(iiiUj)/dxidxj  =  div  (wAu)  +  V2  (u2/ 2) 
to  express  the  solution  (2.1.8)  of  Lighthill’s  equation  in  the  form 
p(x,t)  =  p1(x,t)+p2(x,t), 


(2.3.3) 


Pi  CM)  = 


-p0Xj  d 

47tc0|x|2  dt 


/>Au)(py,  P2(x,()  =  sgR|j/[iU2]d3y,  M-+~. 

(2.3.4) 
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where  [  ]  denotes  evaluation  at  time  t-  |x  -  y|/c0  «  t  -  |x|/c0  +  x.y/c0|x|.  When  retarded 
time  variations  across  the  source  region  are  neglected  the  identity  (2.3.3)  and  the  divergence 
theorem  imply  that  /[w  Au]d3y  =  0,  since  u  ~  0(l/|y|3)  as  |y|  ->•  oo.  To  estimate  the  value 
of  the  first  integral  in  (2.3.4)  it  is  therefore  necessary  to  expand  (u>  A  u)(t  -  |x  -  y|/c0)  in 
powers  of  the  retarded  time  element  x.y/c0|x|.  The  first  term  in  the  expansion  yields, 

p,(x,i)  «  / yi{u>  A  uWi  -  M/O  ~  OW|x|)p„u2M2  (2.3.5) 

where  d/dt  ~  u/l. 

The  order  of  magnitude  of  p2(x,t)  is  estimated  by  first  writing  the  momentum  equation 
in  the  form, 

dn/dt  +  V  (J  dp/p  +  v2/2  +  dcp/dtj  =  —cj  Au-wA  V<p. 

Take  the  scalar  product  with  u 

d(\u2)/dt  +  div  (u (/  dp/p  +  |u2  +  dip/dt)^j  =  — u.w  A  V< p, 
and  integrate  over  the  whole  of  space.  The  contribution  from  the  divergence  vanishes  because 
u(f  dp/p  +  \v2  +  dip/dt )  tends  to  zero  at  least  as  fast  as  l/|y|3  as  |y|  oo.  Hence,  using 
(2.3.2) 

§~t  /  =  ~  I ^U‘W  A  VvP^y’  ^rf3y  ~  ^u3m2>  (2-3-6) 

w'here  the  final  estimate  is  really  a  crude  upper  bound,  that  takes  no  account  of  the  details 
of  the  interactions  between  the  vorticity  and  irrotational  velocity. 

Thus, 

p2(x,t)  ~  O  ((e/\x\)p0u2MA), 

and  by  comparison  with  (2.3.5)  we  see  that  p2  ~  0{M2)p\  pi  when  M  <  1.  The 
component  div(p0w  Av)  of  d2(p0ViVj)/dXidxj  is  therefore  the  principal  source  of  sound  at  low 
Mach  numbers  [10,  11].  This  is  consistent  with  Lighthill’s  hypothesis  that  the  back-reaction 
of  the  sound  on  the  flow  is  negligible,  because  (2.3.6)  states  that  the  rate  of  change  of  the 
kinetic  energy  of  the  source  flow  due  to  compressible  effects  is  at  most  of  order  p0uzM2/i 
per  unit  volume,  which  is  much  smaller  than  the  power  p0uz/i.  needed  to  maintain  the  flow. 

Acoustic  analogy  in  terms  of  the  total  enthalpy.  This  analysis  and  the  discussion  at  the 
end  of  §2.1  show  that  the  dominant  aeroacoustic  sources  at  low  Mach  numbers  are  vorticity 
and  entropy  fluctuations.  In  order  to  cast  Lighthill’s  equation  in  a  form  that  emphasizes 
vorticity  and  entropy  as  the  sound  sources,  it  is  necessary  to  select  the  total  enthalpy  B  = 
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w+^v2  an  the  independent  acoustic  variable  ( w  being  the  specific  enthalpy).  This  reduces  to 
-dipjdt  in  irrotational  regions,  where  is  the  velocity  potential  of  the  unsteady  motion.  B  is 
therefore  constant  in  steady  irrotational  flow,  and  far  from  the  acoustic  sources  perturbations 
in  B  =  -dp/dt  represent  acoustic  waves. 

To  reformulate  Lighthill’s  equation  in  terms  of  B  for  a  homogeneous  fluid,  whose  chemical 
composition  is  the  same  everywhere,  we  start  from  Crocco’s  form  of  the  momentum  equation 

[9] 

dv/dt  +  VB  =  -wAv  +  TVs  +  tr  (2.3.7) 

Oi  =  (1  /p)daij/dxj,  (2.3.8) 


and  re-arrange  the  continuity  equation  (2.1.4)  to  read 


1  Dp  PTDs 

— o*rv7  +  dlv  v  = TT7) 

pC 2  Dt  Cp  Dt 

where  Dp/Dt  has  been  eliminated  by  means  of  the  thermodynamic  relations 


,  dP  . 

dp  -  + 


(dp) 

[dsj 


ds  “d  (I)  =(lf 

V  \  /  P  \ 


),(£), 


_-bpT 


(2.3.9) 


P  and  cp  being  respectively  the  coefficient  of  expansion  and  the  specific  heat  at  constant 
pressure.  For  an  ideal  gas  (3  =  1/T. 

Subtract  the  divergence  of  (2.3.7)  from  the  time  derivative  of  (2.3.9): 


|-  -  V2S  =  div(w  A  v  -  TVs  -  <r)  +  ^  (—j£)  ■  (2-3.10) 

dt  \pc 2  Dt)  v  ’  dt\cp  Dt) 

This  equation  is  already  in  a  form  that  strongly  suggests  that  the  terms  on  the  right  hand  side 
constitute  the  most  important  sources  of  flow  generated  sound,  because  they  are  unchanged 
when  the  compressibility  K4  =  1/pc2  — y  0,  and  accordingly  dominate  the  hydrodynamic  far 
field  of  B  which  must  ultimately  match  the  outgoing  sound  waves. 


The  first  term  on  the  left  of  (2.3.10)  is  expanded  as  follows: 

8_  (j_Dp\  _  £_  (J_dp\  1  dvj  dp  (  d  (  1  \  dp  d_  /j_\  dp\ 

dt  l  pc2  Dt )  Dt  l  pc2  dt )  pc 2  dt  dxj  ^  1  dt  \pc2 )  dxj  dxj  \pc2 )  dt  J 

(2.3.11) 

The  last  term  on  the  right  will  be  discarded.  It  vanishes  when  the  adiabatic  compressibility 
Ks  =  1/pc2  is  a  function  of  the  pressure  alone,  i.e.,  when  the  motion  is  homentropic,  or  in 
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the  particular  case  of  an  ideal  gas  where  Ks  =  I/7P  and  7  =  cp/cv  =  constant.  Variations  of 
7  in  a  real  gas  are  significant  only  in  the  presence  of  large  temperature  gradients;  problems 
of  this  kind  are  best  treated  separately.  In  liquids  the  compressibility  may  be  regarded  as 
constant;  acoustic  sources  involving  variations  of  K4  are  important  only  when  gas  bubbles 
are  present.  The  other  terms  on  the  right  of  (2.3.11)  are  transformed  as  follows.  First: 


1  dp 
p  dt 


1  Dp 
p  Dt 


1  _ 

-v.Vp 

P 


db_ds 
Dt  Dt 


-V.<7, 


(2.3.12) 


where  the  thermodynamic  identity  dw  =  Tds  +  dp/p  and  the  momentum  equation  have  been 
used,  and  second,  from  (2.3.7), 


BHi  BE.  =  _ Vp.  (VS  +  w  A  v  -  TVs  -  a)  .  (2.3.13) 

dt  dxj 


Thus,  substituting  from  (2.3.11)  into  (2.3.10),  we  arrive  at  the  acoustic  analogy  equation 
for  the  total  enthalpy 


Vorticity  and  entropy  sources.  The  operator  on  the  left  of  this  equation  is  nonlinear.  All 
terms  on  the  right  hand  side  vanish  in  irrotational  regions,  and  in  the  absence  of  such  terms 
and  of  boundary  motions,  B  =  constant.  The  radiation  condition  ensures  that  the  terms 
on  the  right  may  be  identified  as  acoustic  sources.  They  are  confined  to  the  region  in  which 
uj  zfi.  0  and  where  Vs  7^  0.  The  wave  operator  describes  propagation  of  the  sound  through 
the  non-uniform  flow  although,  as  for  Lighthill’s  equation,  it  will  not  usually  be  permissible 
to  neglect  the  interaction  of  the  aero  dynamically  generated  sound  with  the  vorticity  and 
entropy  gradients  when  the  source  flow  is  very  extensive.  The  following  special  cases  should 
be  noted: 

High  Reynolds  number,  homentropic  flow:  When  dissipation  is  neglected  and  s  = 
constant  (2.3.14)  becomes 

(i  (?£)  -  ?v-H B  -  ^ A  v)-  (2-3i5) 

At  the  low  Mach  numbers  relevant  in  underwater  applications,  when  the  flow  may  be  assumed 
to  be  at  rest  at  infinity  (where  p  =  p0  and  c  =  c0)  further  simplification  is  possible  by  (a) 
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neglecting  nonlinear  effects  of  propagation  and  the  scattering  of  sound  by  the  vorticity,  and 
(b)  taking  c  =  c0,  and  p  =  p0-  Then 

-  V2)  £  =  div(w /\  v),  (2.3.16) 

and  in  the  far  field  the  acoustic  pressure  is  given  by  the  linearized  approximation 

p(x,  t)  »  p0B(x,  t).  (2.3.17) 

Non-homentropic  source  flow:  When  the  source  region  is  not  homentropic,  it  is  important  to  distinguish 
between  terms  on  both  sides  of  (2.3.14)  that  account  principally  for  scattering  of  the  sound,  and  those  (on 
the  right)  that  can  unequivocally  be  recognized  as  sources.  For  compact  sources  we  can  argue  that  the 
correct  source  terms  are  those  remaining  when  the  fluid  is  temporarily  taken  to  be  incompressible.  These 
are  the  sources  on  the  right  of  (2.3.10).  When  dissipation  is  ignored,  (2.3.14)  then  reduces  to 

{(l  +  U'V)  [?(^  +  UV)]  -iv.(PV)}B  =  div(»Av-7'V*)  +  |(^g).  (2.3.18) 

This  approximation  is  useful  when  the  background  mean  flow  may  be  regarded  as  irrotational  (steady  flow 
past  a  streamlined  body,  for  example)  at  mean  velocity  U(x)  and  density  and  sound  speed  p(x)  and  c(x). 
At  very  small  Mach  numbers  it  simplifies  further  to 

(*£  - 4 B ■  a  v  - TVs) + 1-,  (f  k)  ■  <2-319) 

The  final  entropy  source  on  the  right  of  (2.3.18)  and  (2.3.19)  represents  the  production  of  sound  by 
unsteady  heating  of  the  fluid.  According  to  the  equation  of  continuity  (2.3.9)  this  source  is  equivalent  to  a 
volume  monopole  of  strength, 

<2'3'20) 


§2.4.  Compact  Green’s  Function 

The  compact  Green’s  function  provides  a  formal  and  intuitive  procedure  for  calculating  the 
leading  order  monopole  and  dipole  terms  in  the  multipole  expansion  of  the  sound  produced 
by  sources  near  a  solid  body. 

Green’s  function  is  the  causal  solution  of  the  wave  equation 

(^^2  -  V2)  G(x,y,t-r)  =5(x-y)5(t-r),  G{x,y,t  -  r)  =  0  for  t  <  r;  (2.4.1) 
the  time-harmonic  ( frequency  domain )  Green’s  function  satisfies 


(V2  +  K„)G(x,y;o;)  =  5(x-y),  k0  =  u/c0. 


(2.4.2) 
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For  the  time  harmonic  case  all  quantities  vary  in  proportion  to  e  Mt,  and  G(x,  y;  u)  has 
outgoing  wave  behavior.  The  frequency  and  time  domain  Green’s  functions  are  related  by 

G(x,  y ,  t  -  t)  =  -=i  r  G(x,  y;  w)c-“<‘-'>  du,  (2.4.3) 

I'K  J-oo 

where  the  path  of  integration  in  the  w-plane  runs  above  all  singularities  of  G(x,y;o»). 

The  time-harmonic  Green’s  function  satisfies  the  reciprocal  property  [12] 

G(x,  y;  u)  =  G(  y,  x;  w).  (2.4.4) 


Time-harmonic  problems.  Consider  the  Green’s  function  G(x,  y;  u>)  in  the  presence  of 
an  acoustically  compact  body  of  diameter  t  in  the  particular  case  in  which  the  source  point 
y  is  close  to  S  and  the  observer  at  x  is  in  the  acoustic  far  field.  Let  G  have  vanishing  normal 
derivative  on  S.  The  compactness  condition  implies  that  K0t  <C  1. 

Formula  (2.4.4)  enables  the  problem  of  determining  G(x,y;u>)  to  be  posed  as  a  reciprocal 
scattering  problem  in  which  the  spherical  wave 

_pt/co|x-y| 

G'(y'Xiw)  =  SF^i  (2-4-5) 

generated  by  a  point  source  <5(x  -  y)  at  the  far  field  point  x  is  incident  on  S,  and  G(x,  y;  w)  is 
calculated  by  taking  the  independent  variable  in  (2.4.2)  to  be  y  rather  than  x.  Then,  if  the 
coordinate  origin  is  within  S  and  y  is  close  to  S,  the  condition  K0i  1  permits  G(x,y,uj) 
to  be  expanded  in  the  form  [13  -  15] 


G(x,y,u) 


y  ~  0(£),  |x|  -»  oo.  (2.4.6) 

The  first  term  in  the  large  braces  represents  the  incident  wave  (2.4.5)  evaluated  at  y  =  0. 
The  next  term  is  0(ko£)  and  includes  a  component  —in0Xiyi/\x\  from  the  incident  wave  plus 
a  correction  -z/c0£i<£i(y)/lxl  due  to  S.  To  this  order  of  approximation  V^(y)  =  yi~  <p*( y)  is 
a  solution  of  Laplace’s  equation  satisfying  dYi/dyn  =  0  on  S.  Since  (y)  must  decay  with 
distance  from  S,  this  implies  that  ip*  is  simply  the  velocity  potential  of  the  incompressible 
motion  that  would  be  produced  by  translational  motion  of  S  as  a  rigid  body  at  unit  speed 
in  the  i-direction.  The  remaining  terms  in  (2.4.6)  are  of  order  (k 0£)2  or  smaller.  When  they 
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are  neglected  the  resulting  approximation  for  G  can  be  used  to  determine  the  monopole 
and  dipole  terms  in  the  multipole  expansion  of  the  solution  of  a  fluid-structure  interaction 
problem. 

The  potentials  <p*  are  uniquely  defined  by  the  shape  of  the  body  and  satisfy 

M(y)  =  n,  on  S,  (2.4.7) 

dyn 

where  n  is  the  unit  normal  on  S  directed  into  the  fluid.  When  the  body  is  rigid  they  also 
determine  the  added  mass  tensor  My  [9],  which  is  symmetric  and  given  by 

My  =  -p0  j>s  <P* (y)njdS(y)  =  -p0  ^  (p*d<p*/dyndS(y)  =  -p0  £  ip*d(p*  /  dyndS(y)  (2.4.8) 

A  body  in  translational  accelerated  motion  at  velocity  U(t)  in  an  ideal,  incompressible  fluid 
experiences  a  reaction  force  from  the  fluid  equal  to  —d(MijUj)/dt  in  the  i-direction  [9].  An 
external  force  F{  acting  through  the  center  of  mass  of  a  body  of  mass  m  therefore  produces 
accelerated  motion  determined  by  d(mUi  +  Mjj Uj)/dt  =  Fi. 

Time-domain  problems.  The  general  solution  of  the  inhomogeneous  wave  equation 

Qr  -  v2)  p(x>  f)  =  ^(x-  *)>  (2-4-9) 

in  the  region  outside  a  stationary  surface  S  can  be  written 

p{x,t)  =  ^p{y,r)^-(x,y,t-T)  +  p0^(y,T)G(x,y,t-T)^njdS(y)dT 

+  J  G(x,  y,  t  -  t)F( y,  r)d?ydT,  (2.4.10) 

When  S  is  rigid  we  choose  Green’s  function  to  have  vanishing  normal  derivative  on  S.  Then, 
using  the  expansion  (2.4.5)  to  calculate  G(x,y,t  -  r)  from  (2.4.3),  we  find 

P(x,i)  «  ^  (  / ^(y,  t  -  |x|/c0)d3y  +  /  ^(y)^(y,  t  -  \x\/c0)d3y 

+  E  /  $n(x/lxl>  y)^(y> t  -  lxl/co)<*3y)  >  lxl  ->  °°-  (2.4.11) 

n>2  co  01  J  ' 

The  term  of  order  n  in  the  series  ~  0((u£/c0)n).  When  the  body  is  compact,  so  that 
u>£/c0  -C  1,  it  is  usually  sufficient  to  retain  only  the  first  two  terms  in  the  large  brackets,  the 
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monopole  and  dipole.  This  is  equivalent  to  approximating  Green’s  function  by 

1  (  x.Y  \ 

G(x,  y,  t  -  t)  «  -  r  -  |x|/c0)  +  -r-  |x|/c0)  J 

~  An\x\5^~T  ~  dXl  _X-Y/IXD/C°) 

w  - - -S(t-T-  |x-  Y|/c0),  |x|->oo.  (2.4.12) 

47T  x  -  Y|  v  ' 


This  result  can  be  made  symmetric,  in  accordance  with  reciprocity,  by  replacing  x  by 
X  =  x  —  v?*(x),  after  which  we  define  the 

Compact  Green’s  function  for  a  body  bounded  by  a  surface  S: 

G(x,y,t “  T)  =  felX-YI^  " T  "  |X “  Y|/c»)’  (2-4’13) 

where  X{  -  Xi  -  (x),  Yi  =  yi~  (p*{  (y) ,  and  ip*  is  the  velocity  potential  of  the  incompressible 

flow  that  would  be  produced  by  rigid  body  motion  of  S  at  unit  speed  in  the  i-direction,  so 
that  dG(x,  y,t-  T)/dxn  =  dG(x,  y,  t  -  r)/dyn  =  0  on  S. 

In  the  frequency  domain  the  corresponding  approximation  is 


G(x,y;w)  = 


_eiKo  |X-Y| 

4tt|X  -  Y] 


(2.4.14) 


Note  that,  because  of  the  symmetrical  manner  in  which  x  and  y  occur  in  (2.4.13)  and 
(2.4.14),  the  location  of  the  coordinate  origin  is  now  arbitrary.  However,  it  is  still  necessary 
that  either  x  or  y  should  lie  in  the  far  field  of  the  body.  When  both  x  and  y  are  in  the 
far  field  predictions  made  with  the  compact  Green’s  function  will  be  the  same  as  when  the 
body  is  absent.  This  is  because,  for  distant  sources  the  amplitude  of  the  sound  scattered 
by  a  compact  rigid  object  is  0((ko£)2)  smaller  than  the  incident  sound,  i.e.,  is  of  quadrupole 
intensity  [12]. 
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§2.5.  Surface  force  exerted  by  a  rigid  body 

Force  exerted  on  an  incompressible  fluid.  At  low  Mach  numbers  the  sound  generated 
by  turbulent  flow  interacting  with  a  compact  body  depends  on  the  surface  forces  calculated 
as  if  the  fluid  is  incompressible  (see  §2.6).  The  force  F(t)  exerted  on  the  fluid  by  a  solid  in 
rigid  body  motion  at  velocity 

U  =  U0  +  0  A  (x  —  x0{t)) , 

where  U0  =  dx0/dt  is  the  velocity  of  its  center  of  volume  x0(t),  and  fl(t)  is  its  angular 
velocity  (Figure  2.1),  is  given  by  the  following  general  formula  in  terms  of  U0  and  w  [16] 

F(()41/vftxWx_m',ir’  (251) 

where  m0  is  the  mass  of  the  fluid  displaced  by  the  body.  This  equation  is  applicable  for 
three  dimensional  bodies.  In  two  dimensions  m0  is  the  mass  displaced  per  unit  span  and  the 
factor  of  |  is  omitted. 

The  integral  in  (2.5.1)  defines  the  impulse  of  the  coupled  system  [9],  which  is  an  invariant 
of  the  motion  when  the  body  is  absent.  To  evaluate  the  integral  the  vorticity  must  be 
defined  .everywhere  including  the  region  occupied  by  the  body,  i.e.,  vortex  lines  must  be 
continued  into  the  interior  of  S  to  form  re-entrant  filaments.  For  a  non-rotating  body,  vortex 
lines  meeting  the  surface  form  re-entrant  loops  by  continuation  on  the  surface.  The  impulse 
then  becomes  a  function  of  time  because  the  motion  of  these  vortex  lines  within  and  on  S 
is  no  longer  governed  by  the  Navier-Stokes  equations.  For  irrotational  flow  u>  =  0  in  the 
fluid,  but  not  on  S,  where  it  must  be  taken  as  the  singular  distribution  of  bound  vorticity 
Wj  =  nA  (V<p  -  U)5(f),  where  f(x,t)  =  0  defines  the  surface  S  moving  at  velocity  U,  and 
<p  is  the  velocity  potential  of  the  irrotational  flow;  inside  S  we  must  take  oo  —  20. 

At  high  Reynolds  numbers,  involving  the  interaction  of  turbulent  flows  with  S,  the  surface 
boundary  layers  are  often  very  thin,  and  the  exterior  fluid  motion  can  frequently  be  approxi¬ 
mated  by  irrotational  flow  with  a  superposed  distribution  of  vorticity.  In  these  circumstances 
it  would  be  convenient  to  have  an  expression  for  F  that  minimizes  the  contribution  from  the 
bound  vorticity.  Such  a  formula  can  be  derived  for  a  body  in  translational  motion  [17  -  19]. 

To  do  this,  suppose  the  fluid  is  at  rest  at  infinity,  and  write  the  force  in  the  form 

Fi  =  4:  [  p0Vid3x  +  <f  pnidS,  (2.5.2) 

dt  Jv  J s 
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E 


Figure  2.1.  Motion  of  a  rigid  body  in  fluid  at  rest  at  infinity. 
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where  V  here  denotes  the  fluid  between  S  and  a  large  fixed  control  surface  E  (Figure  2.1) 
whose  normal  n  is  directed  into  V.  To  express  the  right  hand  side  of  this  formula  in  terms  of 
the  vorticity  and  the  translational  velocity  of  the  body,  we  first  cast  the  momentum  equation 
for  incompressible  flow  in  the  form 

+  V  (p  +  =  -Po“>  A  v  -  rjcuvl  w.  (2.5.3) 

The  scalar  product  of  this  equation  with  V-X^x,*)  =  V  (a*  -  <p-(x,t))  is  then  integrated 
over  V.  Xi  is  the  function  used  in  §2.4  to  define  the  compact  Green’s  function;  it  is  here 
dependent  on  time  because  of  the  translation  of  S.  By  applying  the  divergence  theorem,  the 
integrated  equation  can  be  manipulated  to  yield  an  alternative  representation  of  the  right 
hand  side  of  (2.5.1).  Using  the  definition  (2.4.8)  of  the  added  mass  tensor,  the  following 
formula  is  then  obtained  for  the 

Force  exerted  on  an  incompressible  fluid  by  a  body  in  translational  motion 

Fi  =  Mi,—  -  po  [  VXj.w  A  Vre;d3x  +  T]  <f  u  A  VXj.dS,  (2.5.4) 

J  dt  Jv  J  s 

where  vre/  =  v  -  U  is  the  fluid  velocity  relative  to  S.  Bound  vorticity  makes  no  contribution 
to  the  volume  integral  because  (i)  vre/  =  0  on  S  for  viscous  flow,  or  (ii)  VIi.wAvre,  =  0  on 
S  for  an  ideal  fluid,  since  VXit  u>  and  vre(  are  then  all  locally  parallel  to  the  surface. 

This  formula  expresses  the  force  as  a  sum  of  three  essentially  distinct  components:  (i)  the 
inertia  due  to  the  added  mass  of  the  body,  (ii)  the  vector  sum  of  the  normal  stresses  induced 
on  S  by  the  vorticity,  (iii)  the  viscous  skin  friction. 

Force  exerted  by  a  semi-infinite  surface.  The  unsteady  normal  force  F2  produced  by 
turbulence  near  an  infinitely  long,  stationary,  rigid  strip  (Figure  2.2a)  of  chord  2a  is  given 
by  (2.5.4)  by  taking 

X2  =  Re(  -  iVz2  -  a2),  z  =  xi+ix2 ,  (2.5.5) 

where  the  strip  occupies  |a;i|  <  a,  x2  =  0,  —00  <  X3  <  00. 

Consider  vorticity  confined  to  the  neighborhood  of  the  edge  X\  =  +a  of  the  strip  within 
a  distance  ~  t  from  the  edge.  When  t  <  2a,  X2  can  be  approximated  within  the  vortical 
region  by 

X2  «  V2a  <p*(x),  (2.5.6) 


where  (Figure  2.2b) 


<p*(x)  =  v^sin  ,  X\—  a  =  rcos0,  X2  =  rsinP.  (2.5.7) 
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Figure  2.2.  Force  generated  by  vortex-edge  interaction. 
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y*(x)  is  just  the  velocity  potential  of  an  incompressible  flow  around  the  edge  (xx  =  +a)  of  a 
semi-infinite  rigid  half  plane  extending  to  xx  =  -oo.  The  order  of  magnitude  of  the  surface 

force  is  then  _ 

F2  =  -V2 a  I  p0(u  A  v)(x,  t).V<p*(x)d3x  ~  p0v2i2\jj  , 

which  is  unbounded  as  a/ 1  — >•  oo. 

But  in  this  limit  the  airfoil  chord  ultimately  becomes  non-compact,  and  it  is  necessary 
to  take  account  of  compressibility  in  calculating  the  surface  force,  in  which  case  it  might 
be  anticipated  that  the  role  of  the  chord  length  2 a  in  determining  the  surface  force  is  then 
played  by  the  characteristic  acoustic  wavelength  of  the  aerodynamically  generated  sound.  If 
the  edge  region  occupied  by  the  vorticity  is  acoustically  compact,  the  force  can  be  calculated 
by  making  use  of  the  approximation  given  below  in  equation  (4.1.4)  to  Green’s  function  for 
a  rigid  halfplane,  which  yields  the: 


Frequency  domain  force 

t7f  /  .  17T 

F2{u)  «  -^==Jpo(w  Av)(x,w).ViJ*(x)rf3x  =  — ft(wAv)(x,w).V^(x)rf3x, 
^VK°  (2.5.8) 

where  A  =  2i t/k0  is  the  acoustic  wavelength; 

Time  domain  force 


l^/Co 
7T 


LdrI 


(oj  A  v)(x,t).V</?* (x) 
\/i  —  r 


orx. 


(2.5.9) 


These  formulae  remain  valid  when  the  edge  of  the  half  plane  supports  a  rigid,  “attachment” 
(Figure  2.2c)  along  its  span  whose  cross-section  is  acoustically  compact  and  whose  length 
scale  of  variation  in  the  spanwise  (£3)  direction  is  compact.  In  this  case  the  potential  function 
< p*(x )  =  (p*(xi,x2)  must  be  replaced  by  <h*(x)  =  $*(:ri,  £2,  £3)  which  describes  potential  flow 
around  the  irregular  edge,  and  satisfies 

$*(£!,  £2,  x3)  “*■  ^*(^1.^2)  =  Vrsin  ^  ,  as  r  =  (x[  -  a)2  +  ~x\  ->•  00. 
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§2.6  Sound  generation  in  the  presence  of  rigid  boundaries 
The  Ffowcs  Williams  -  Hawkings  equation  (2.2.4)  enables  the  aerodynamic  sound  to  be 
expressed  as  the  sum  of  the  sound  generated  by  quadrupole  sources  together  with  contribu¬ 
tions  from  monopole  and  dipole  sources  distributed  on  boundaries.  For  turbulent  flow  near 
a  fixed  rigid  surface,  the  direct  sound  from  the  quadrupoles  Tij  is  augmented  by  radiation 
from  surface  dipoles  whose  strength  is  the  force  per  unit  surface  area  exerted  on  the  fluid.  If 
the  surface  is  in  accelerated  motion  there  are  additional  dipoles  and  quadrupoles,  and  neigh¬ 
boring  surfaces  in  relative  motion  produce  “potential  flow”  interactions  that  generate  sound. 
At  low  Mach  numbers  the  acoustic  efficiency  of  the  surface  dipoles  exceeds  the  efficiency  of 
the  volume  quadrupoles  by  a  large  factor  ~  0(1/M2).  Thus,  the  presence  of  solid  surfaces 
within  low  Mach  number  turbulence  can  lead  to  substantial  increases  in  aerodynamic  sound 
levels. 


Acoustically  compact  bodies.  Consider  the  production  of  sound  by  turbulence  near  a 
compact,  stationary  rigid  body.  Let  the  fluid  have  uniform  mean  density  p0  and  sound  speed 
c0,  and  assume  the  Mach  number  is  sufficiently  small  that  convection  of  the  sound  by  the 
flow  may  be  neglected.  The  acoustic  pressure  in  the  far  field  is  given  by  Curie’s  equation 
(2.2.5)  in  the  form 


p(x,  t) 


XiX 


1^3 


d2 


47tc2|x|3  dt2 


I  [poViV^y  +  js  bp«]<*S(y),  |x|  ->  oo,  (2.6.1) 


where  the  terms  in  square  brackets  are  evaluated  at  the  retarded  time  t  —  |x  —  y|/c0,  p[j  = 
{p  ~  Po)$ij  —  &ij  is  the  compressive  stress  tensor  and  the  unit  normal  n  on  S  is  directed  into 
the  fluid.  The  first  term  on  the  right  is  the  quadrupole  noise,  which  has  the  same  formal 
structure  as  when  the  body  is  absent.  The  surface  integral  is  the  dipole  contribution. 


In  practice  the  length  scale  i  of  the  turbulence  is  comparable  to  some  dimension  of  the 
body.  The  frequency  of  the  hydrodynamic  fluctuations  ~  v/ft,  and  the  whole  motion  in  the 
neighborhood  of  the  body  is  therefore  acoustically  compact  if  ( v/l)Zfc0  =  v/cQ  <C  1,  i.e., 
provided  the  Mach  number  is  small.  The  retarded  time  variations  on  S,  which  are  of  order 
x.y/c0|x|  when  |x|  — >  oo  may  then  be  neglected  so  that,  with  the  origin  in  the  neighborhood 
of  the  body,  the  first  approximation  to  the  surface  integral  yields  the  dipole  sound  field 


Xj  dFj 
47tc0|x|2  dt 


(t  -  |x|/c0), 


where  Fi(t)  =  fsp'j(y,  t)rijdS(y)  is  the  net  force  exerted  on  the  fluid  by  the  rigid  body. 
In  this  compact  approximation  the  force  can  be  calculated  in  terms  of  the  vorticity  u>  and 
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the  velocity  v  from  equation  (2.5.4)  (in  which  dUj/dt  =  0  for  a  stationary  body),  and  the 
acoustic  pressure  cast  in  the  form 

p(x,f)  «  ( J  [(w  A  v).  VF,]  <f3y  ~  ^  (u)  A  VF.n]  rfS(y)} 

S3  j^|j|/(«Av)(y,(-|x|/c0).vyi(y)rf3y.  M -><*>,  (2.6.2) 

where  [  ]  denotes  evaluation  at  t  -  |x|/c0,  and  Yt  =  yt  -  <p*(y)  is  the  velocity  potential  of 
incompressible  flow  past  the  body  having  unit  speed  in  the  i-direction  at  large  distances  from 
S  (the  corresponding  contribution  to  the  radiation  being  equivalent  to  a  dipole  orientated 
in  the  i-direction) .  The  approximation  in  the  second  line  is  applicable  at  high  Reynolds 
numbers,  vi/u  >  1,  when  the  “skin  friction”  contribution  to  the  surface  force  is  negligible. 

The  order  of  magnitude  of  F  ~  p0v2£2,  and  d/dt  ~  v/L  The  dipole  generated  acoustic 
pressure  is  therefore  of  order  (£/|x| )p0v2M  ( M  =  v/c0),  with  corresponding  sound  power 
na  ~  £2p0v3M3  and  acoustic  efficiency  Yla/(£2p0v3)  0(M3).  For  the  direct  quadrupole 

radiation,  the  acoustic  pressure  ~  (£/|x| )pav2M2,  the  same  as  in  the  absence  of  the  body, 
and  the  efficiency  is  ~  0(M5).  The  radiation  is  accordingly  dominated  by  the  dipole  when 
M  is  small,  and  as  M  — >  0  the  acoustic  power  exceeds  the  quadrupole  power  by  a  factor 
~  1/M2  »  1.  Precisely  how  small  M  should  be  for  this  to  be  true  depends  on  the  details  of 
the  flow. 

This  increase  in  acoustic  efficiency  brought  about  by  surface  dipoles  on  an  acoustically 
compact  body  occurs  also  for  arbitrary,  non-compact  bodies  when  vorticity  interacts  with 
compact  structural  elements,  such  as  edges,  corners  and  protuberances.  In  the  extreme  case 
of  a  flat ,  rigid  surface,  whose  dimension  greatly  exceeds  the  acoustic  wavelength,  the  normal 
force  exerted  on  an  incompressible  fluid  vanishes  identically.  The  integrated  contribution 
from  the  dipoles  is  then  null  unless  account  is  taken  of  retarded  time  differences,  when 
their  aggregate  effect  becomes  the  same  as  a  higher  order  quadrupole  that  coincides  with 
the  image  of  the  turbulence  volume  quadrupoles  in  the  wall  [20];  the  overall  efficiency  is 
therefore  0(M5),  the  same  as  for  quadrupoles  in  free  space. 

Application  of  compact  Green’s  functions.  These  results  can  also  be  obtained  and 
extended  by  making  use  of  the  compact  Green’s  function  of  §2.4.  Consider  vorticity  convected 
past  the  body  in  low  Mach  number,  homentropic  mean  flow  (Figure  2.3). 
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Figure  2.3.  Turbulent  flow  past  a  stationary,  rigid  body. 
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The  acoustic  radiation  will  be  determined  by  solving  equation  (2.3.16).  The  outgoing  wave 
solution  is  given  by  Kirchhoff’s  formula  [21] 

B(x,t)  =  £|B(y,T)|^(x,y,t-r)-|^(y,r)G(x,y,t-r)|nidS(y)dr 

+  J  G(x,  y,  t  -  r)(div(u;  A  v)(y,  r))d3ydr,  (2.6.3) 

where  the  retarded  time  integration  with  respect  to  r  is  taken  over  (-oo,  oo). 

The  Green’s  function  is  chosen  to  have  vanishing  normal  derivative  on  S,  and  the  volume 
integral  is  transformed  by  the  divergence  theorem: 

f  Gdiv(u>  A  v)d3y  =  —  j>s  G(oj  Av)jUjdS  —  f(u>  Av).VGd3y. 

When  viscous  effects  on  S  are  assumed  to  be  the  same  as  for  an  incompressible  fluid,  Crocco’s 

equation  (2.3.7)  may  be  approximated  by 

dv/dt  +  VB  —  — oj  A  v  -  v  curl  u> 

and  (2.6.3)  becomes 

B(x,t )  «  -J  (w  Av)j(y,r)-  (x,y,t  -  r)d3ydr 

£  w(y,  r)  A  |(x,  y,  t  -  r).dS{y)dT  +  jf  G(x,  y,  t  -  r)^( y,  r)njdS(y)dr 
7  (2.6.4) 

In  the  absence  of  mean  flow  B  =  p/p0  in  the  acoustic  far  field.  This  remains  true  in  the 
presence  of  flow  at  very  low  Mach  number  M0,  since  in  general  the  pressure  as  |x|  oo  is 
calculated  from 

%~Po(]ti  +  U°-V)  B  =►  P(x,*)«PoS(x,t)(l  +  0(M0)).  (2.6.5) 

The  last  term  on  the  right  of  (2.6.4)  vanishes  when  S  is  rigid.  Thus,  when  the  G  is 
approximated  by  the  compact  Green’s  function  (2.4.13)  we  find,  correct  to  dipole  order, 

-  /[(«AV)-^.  (2-6'6) 
where  [  ]  denotes  evaluation  at  t  —  |x|/c0;  in  the  second  line  the  skin  friction  contribution 
has  been  discarded  for  high  Reynolds  number  flow.  These  formulae  are  equivalent  to  Curie’s 
approximations  (2.6.2). 
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If  the  body  also  executes  small  amplitude  vibrations  at  velocity  v(t),  the  compact  Green’s 
function  used  in  (2.6.4)  yields 


p(x,t) 


Xi 


47tc0|x|2  dt 


“b  ttIq 


dvi 

dt 


,  |x|  — >■  oo, 


(2.6.7) 


J  t-\x\/Co 


where  F  is  the  force  exerted  on  the  fluid  and  m0  is  the  mass  of  fluid  displaced  by  the  body. 


By  writing  v  =  U  +  u  in  (2.6.6),  where  U  is  the  irrotational  velocity  component  of  the 
mean  flow  past  the  body,  the  velocity  u  is  identified  with  the  vorticity  u>  =  curl  u.  For 
streamlined  bodies  it  frequently  happens  that  u<U.  No  sound  is  generated  when  u  =  0 
but,  in  contrast  to  sound  production  by  the  volume  quadrupoles  or  when  U  =  0,  mean 
flow  convection  of  vorticity  past  the  body  increases  the  amplitude  of  the  sound  from  being 
proportional  to  u2  to  uU.  This  linear  component  of  the  radiation  is  furnished  by  the  source 
div(a>AU).  Furthermore,  if  U  -»  U„  at  large  distances  from  the  body,  then  ijVhi.wAU  =  0 
when  x  is  parallel  to  U0  because  U(y)  =  UoiVYi,  so  that  the  linear  surface  dipoles  have  no 
net  component  parallel  to  the  incident  mean  stream  U0. 

For  a  bluff  body  the  magnitude  of  the  rotational  velocity  u  in  the  wake  is  comparable 
to  the  mean  flow  speed,  although  for  moderate  Reynolds  numbers  the  translational  velocity 
of  a  typical  vortical  structure  remains  small  until  it  is  released  into  the  flow.  When  this 
happens  the  core  accelerates  to  a  streamwise  velocity  close  to  that  of  the  mean  flow  and  a 
sound  pulse  is  radiated.  The  periodic  nature  of  these  events  for  wires  in  a  wind  is  responsible 
for  jEolian  tones  [22  -  24].  The  conclusion  above  that  the  intensity  of  the  sound  is  strongest 
in  directions  at  right  angles  to  the  mean  flow,  corresponds  to  dipole  radiation  produced  by 
a  fluctuating  lift  force,  which  reverses  its  direction  as  successive  vortices  are  shed  into  the 
wake. 
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§2.7  Reciprocity  applied  to  sources  adjacent  to  a  compliant  wall 

The  generation  of  sound  by  turbulent  flow  over  an  inhomogeneous  elastic  surface  involves 
the  excitation  of  structural  vibrations  and  their  subsequent  scattering.  The  general  problem 
is  amenable  to  simple  analytical  treatment  only  when  convection  by  the  mean  flow  is  ignored, 
which  is  usually  the  case  in  underwater  acoustics.  The  influence  of  small  amplitude  compli¬ 
ance  can  then  be  estimated  by  the  following  form  of  the  reciprocal  theorem  [12],  which  we 
state  in  a  form  suitable  for  application  to  surfaces  for  which  the  compliant  motion  is  normal 
to  the  undisturbed  surface: 

p(x,  co)gR(uj)  =  -ica(R(y,  x;  u)F(  y,  u)  (2.7.1) 

where  p(x,  u)  is  the  perturbation  pressure  at  a  point  x  in  the  fluid  generated  by  a  normal 
force  F(y,u>)  applied  to  the  surface  (directed  into  the  surface)  at  y,  and  CH(y,x;u>)  is  the 
normal  displacement  of  the  plate  (directed  into  the  fluid)  at  y  produced  by  a  point  volume 
source  of  strength  qR(u>)  at  x. 

This  formula  is  usually  applied  to  the  problems  involving  the  interaction  of  turbulence 
“gusts”  with  an  elastic  surface  in  the  vicinity  of  a  structural  inhomogeneity.  By  taking  the 
volume  source  to  have  unit  strength  and  the  force  to  be  the  gust-generated  blocked  surface 
pressure  ps  (i.e.,  the  surface  pressure  that  would  be  generated  by  the  gust  if  the  surface  is 
regarded  as  rigid),  the  acoustic  pressure  for  each  frequency  u>  is  determined  by  the  equation 

p(x,w)  =  -iw£pi(y,o;)C/i(y,x,a;)dS(y),  (2.7.2) 

where  £R(y ,u)  is  the  flexural  displacement  at  y  produced  by  a  unit  volume  source  at  x.  In 
real  form,  this  becomes 

p(x,  t)  =  ^£^/ooP*(y’rKR(  y,x,t-r)dS(y)dr.  (2.7.3) 
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3.  SOUND  GENERATED  BY  FLOW  OVER  A  FLAT  SURFACE 

Intense  surface  pressure  fluctuations  beneath  a  turbulent  boundary  layer  generate  sound 
and  structural  vibrations.  Sound  is  produced  directly  by  aerodynamic  sources  within  the 
fluid,  and  indirectly  by  the  interaction  of  hydrodynamic  pressures  and  the  surface  vibrations 
with  discontinuities  of  the  wall. 

§3.1  Wall  pressure  wavenumber-frequency  spectrum 

The  blocked-pressure  is  the  pressure  developed  beneath  a  boundary  layer  on  a  hard  wall,  and 
is  twice  the  pressure  that  a  nominally  identical  flow  would  produce  if  the  wall  were  absent. 
The  response  of  a  flexible  wall  to  boundary  layer  forcing  depends  on  both  the  temporal  and 
spatial  characteristics  of  the  pressure.  When  the  wall  is  locally  plane  and  the  fluctuations  can 
be  regarded  as  statistically  stationary  in  time,  the  characteristics  are  expressed  in  terms  of  the 
wall-pressure  wavenumber-frequency  spectrum  P(k,w)  (the  Fourier  transform  of  the  space- 
time  correlation  function  of  the  wall  pressure).  By  convention,  coordinate  axes  (xux2,x3) 
are  taken  with  aq  and  x3  respectively  parallel  and  transverse  to  the  mean  flow,  and  with  x2 
measured  outwards  from  the  wall  (Figure  3.1a).  Then 

P(k,u)  =  75irT  f  k  =  (h,0,h),  (3.1.1) 

[Z7T )°  J -oo 

where 

fc(Vuya,T)  =  {ps(xi,xz,t)ps{xl  +  yux3  +  y3,t  +  r))  (3.1.2) 

is  the  wall  pressure  correlation  function,  ps(x i,  x3,  t)  is  the  fluctuating  component  of  the  wall 
pressure,  and  the  angle  brackets  (  )  denote  an  ensemble  average.  The  wavenumber  k  is 
parallel  to  the  wall,  and  (because  TZ  is  real)  P(-k,  -co)  =  P(k,o>). 

It  is  usually  assumed  that  variations  of  TZ(yi,y3, r)  and  P(k,  u>)  with  x\,  x3  are  significant 
only  over  distances  that  are  much  larger  than  the  length  scales  of  the  dominant  turbulent 
motions,  although  these  distances  may  not  be  large  compared  to  the  wavelengths  of  the 
boundary  layer  generated  sound.  When  TZ(yi,yi,  t)  is  independent  of  (aq,  £3)  and  t,  the  wall 
pressure  is  stationary  random  in  position  and  time,  and 

(ps(k,  w)P;(kV)>  =  S(  k  -  k')<5(u/  -  o/)P(k,u>),  (3.1.3) 

where  the  asterisk  denotes  complex  conjugate. 
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*2 


Figure  3.1.  Turbulent  boundary  layer  and  wall  pressure  spectrum  at 
low  Mach  number  when  cj5/U  1. 
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The  wall  point  pressure  frequency  spectrum  $pp(a>)  is  defined  by 

(p2s{xi,x3,t))  =  r  %p{u)du,  $pp(w)=/  ?(\z,u)dkxdkz.  (3.1.4) 

J— oo  J —oo 

The  root-mean-square  wall  pressure  fluctuation  prms  is  proportional  to  the  dynamic  pres¬ 
sure  \p0U2  {po  and  U  being  the  mean  density  and  velocity  of  the  free  stream) 

Prms  =  =  O  \p0U2\  (3-1.5) 

for  a  smooth,  hard  wall  at  low  subsonic  velocities  a  «  0.01  [25,  26]. 

The  characteristic  shape  of  the  wavenumber-frequency  spectrum  at  low  Mach  numbers  is 
illustrated  in  Figure  3.1b  for  a  fixed  frequency  satisfying  u8/U  >  1,  where  5  is  the  boundary 
layer  thickness  (the  distance  from  the  wall  at  which  the  mean  flow  velocity  is  0.9917).  The 
strongest  pressure  fluctuations  occur  within  the  convective  ridge  centered  on  ki  =  lo/Uc,  kz  = 
0,  where  Uc  is  a  convection  velocity,  and 

P(k,u;)~10  3p2oV6Ju,\  (3.1.6) 

v „  being  the  friction  velocity.  The  velocity  Uc  can  be  interpreted  as  the  translational  speed  of 
the  principal  boundary  layer  eddies,  and  both  experiment  and  numerical  simulations  indicate 
that  Uc  ~  0.5 U  -  0.7 U,  with  only  a  weak  dependence  on  frequency  u. 

At  low  subsonic  Mach  numbers,  the  phase  velocity  ujk  of  surface  pressures  in  the  convec¬ 
tive  domain  is  subsonic  and  these  pressures  therefore  decay  rapidly  with  distance  from  the 
wall,  and  do  not  correspond  to  sound  waves.  However,  convective  pressures  can  generate 
sound  (and  structural  vibrations)  when  the  wall  is  rough  or  has  other  discontinuities  at  which 
convective  energy  is  scattered.  For  a  highly  compliant  elastomeric  wall,  in  which  the  shear 
wave  speed  is  comparable  to  U,  the  unsteady  surface  shear  stresses  that  accompany  the  con¬ 
vective  pressures  may  also  cause  significant  structural  excitation.  For  a  typical  homogeneous 
surface,  the  strongest  coupling  between  the  wall  and  flow  is  usually  attributed  to  longer  wave¬ 
length,  “subconvective”  pressures  in  the  low-wavenumber  region  |«0|  <  k  <C  \u\/Uc  where 
“resonant”  forcing  of  the  wall  occurs  because  the  phase  velocity  u/k  coincides  with  that 
of  a  vibrational  mode  of  the  wall.  Here  the  spectral  levels  are  usually  30  to  60  dB  below 
the  convective  peak.  In  the  acoustic  domain  k  <  |«0|,  the  phase  velocity  is  supersonic,  and 
wall  pressure  fluctuations  are  actually  sound  waves  produced  directly  by  boundary  layer 
quadrupoles.  The  secondary  circular  ridge  in  Figure  3.1b  at  k  =  |/€0|  is  discussed  at  the  end 
of  §3.2. 
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§3.2  Smooth  wall  spectra  at  low  Mach  numbers 

Theoretical  models  of  P(k,  u>)  generally  assume  the  flow  to  be  homentropic,  of  low  Mach 
number  over  a  flat,  rigid  wall  with  no  mean  pressure  gradient.  The  steady  growth  of  the 
boundary  layer  in  the  streamwise  direction  is  usually  ignored,  and  the  mean  flow  velocity  is 
taken  to  be  parallel  to  the  wall  and  dependent  only  on  distance  x2  from  the  wall.  A  sufficient 
body  of  experimental  data  is  available,  derived  principally  from  wind-tunnel  experiments, 
to  permit  the  formulation  of  empirical  models  of  the  convective  domain,  but  reported  prop¬ 
erties  of  the  subconvective  and  acoustic  domains  often  vary  significantly  between  different 
experiments. 

The  hydrodynamic  domain.  Chase  [27]  has  devised  the  following  model  for  P(k,w)  at 
low  Mach  numbers  in  the  hydrodynamic  domain  k  »  |/c0|: 


P(k,qQ 

f%v*8 3 


1 

[(M)2  + 1/&2]1 


(cm(M)2  +  Ct(M)2 


(M)2  + 1  /b2\ 

(k5)2  +  l/b2  ) 


k\  =  (w  -  Ucki)2 / (hvt)2  +  k2, 

M  =  U/c0  <  1,  u5/U>  1.  (3.2.1) 


The  adjustable  coefficients  in  these  expressions  are  fixed  by  comparison  with  experiment, 
and  recommended  values  are 


b  *  0.75,  CM  «  0.1553,  CT  »  0.0047,  h  «  3. 


(3.2.2) 


A  simpler  approximation  which  is  applicable  for  uS/U  >  1  in  the  immediate  neighborhood 
of  the  convective  ridge  is  due  to  Corcos  [28]: 

P(k,u;)  =  +  g2^_  w/[/c) 2] 


i\  «  9Uc/u),  «  lAUc/u), 


(3.2.3) 


where  £1  and  £3  are  turbulence  correlation  lengths  in  the  Xi-  and  ^-directions.  The  point 
pressure  spectrum  $pp,  based  on  data  collated  by  Chase  [27],  may  be  approximated  by 


{PoVlY 


(“S./V)2  „ 


0.12. 


(3.2.4) 
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In  this  formula  <5*  =  J0°°  (1  -  vl{x2)/U)dx2  (~  8/S)  is  the  boundary  layer  displacement 
thickness,  where  Vi(x2)  is  the  mean  velocity  in  the  streamwise  direction  at  distance  x2  from 
the  wall. 

Equations  (3.2.1)  imply  that  P(k,u)  /  p2av353  ~  (kd)2  when  k5  <  1  (the  Kraichnan-Phillips 
Theorem).  In  the  interval  1  <  k5  <  \v\5/Ue  (where  the  term  in  CT  is  still  dominant  in 
(3.2.1))  the  spectrum  is  wavenumber-white,  i.e. 

P(k,u J)/Pyj3  «  C Th3(v*/U)3/{u5/U)3,  l<k6<  \u>\ S/Uc.  (3.2.5) 

This  may  be  contrasted  with  a  low  wavenumber  model  due  to  Sevik  [29],  based  on  measure¬ 
ments  of  P(k,  a>)  at  k  =  0  for  both  smooth  and  rough  walls  on  buoyant  bodies  in  water  and 
in  quiet  wind  tunnels: 

P(0 ,w)/p20vl83  »  127 M2(v*/U)/(lo5/U)4*, 

24  <  \u\5/U  <  240,  0.01  <  M  <  0.15,  M  =  U/c0.  (3.2.6) 

Sevik’s  formula  explicitly  involves  fluid  compressibility  because  of  its  dependence  on  Mach 
number  M,  and  was  originally  proposed  as  a  model  for  the  acoustic  domain.  However,  prac¬ 
tically  all  published  data  suggest  that  the  wall  pressure  spectrum  is  ultimately  wavenumber 
white  as  k5  -»  0,  starting  at  values  of  k  that  are  much  greater  than  the  acoustic  wavenumber 
k0.  If  equation  (3.2.6)  is  applied  at  subconvective  wavenumbers  (1  <  kS  <C  \u}\S/Uc)  pre¬ 
dictions  of  P(0,  oj)  agree  in  magnitude  with  (3.2.5)  at  w5/U  «  24  for  M  =  0.01.  At  higher 
Mach  numbers,  estimates  based  on  Sevik’s  formula  can  be  larger  by  up  to  10  dB,  although 
this  is  offset  by  a  more  rapid  decrease  with  increasing  frequency. 

The  acoustic  domain.  Consider  the  solution  of  Ligh thill’s  equation  expressed  in  Curie’s 
form  (2.2.6)  for  a  boundary  layer  of  infinitesimal  Mach  number,  when  S  coincides  with  the 
plane  of  the  wall.  The  surface  integral  in  (2.2.6)  represents  the  reflection  of  the  quadrupole 
generated  sound  by  the  wall.  When  viscous  shear  stresses  on  the  wall  are  ignored  the 
reflection  coefficient  is  unity,  and  the  pressure  ps  —  c2(p  —  p0)  on  the  wall  x2  =  0  becomes 

ps(x,  t)  =  J  F(y,  t-\x-  y  1  fc0)  Y_  ,  x=(x!, 0,2:3)  (3-2.7) 

where  F(x.,t)  ~  d2 (p0ViVj) / dx,dxj  for  homentropic  flow,  and  the  integration  is  over  the 
unsteady  sources  in  the  boundary  layer. 

When  this  expression  is  used  in  (3.1.1)  to  evaluate  P(k,w),  the  resulting  multiple  integral 
diverges  for  k  =  |/c0|,  at  the  edge  of  the  acoustic  domain,  if  the  boundary  layer  is  assumed 
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to  be  of  unlimited  extent.  The  singularity  is  caused  by  the  collective  action  of  “grazing” 
acoustic  waves,  whose  mean  square  decay  like  l/|x|2  with  propagation  distance  over  the  wall 
is  countered  by  the  growth  in  the  surface  area  of  the  boundary  layer  as  |x|  -»  oo. 

Chase  [27]  has  extrapolated  (3.2.1)  into  the  acoustic  domain,  obtaining 

.P(k,w)  =  1  r  Cuihsye 

pyy  ~  [(M)2+i/62]uifc2-«oi+^ 


+  Cr{k8)2 


{k+S)2  +  1/b2 
(kS)2  +  1/b2 


+ 


czk 2  \ ' 

\k2-K2\  +  e2K2)\’ 


(3.2.8) 


where  cj  =  |,  c2  =  C3  =  |,  and  the  values  of  the  other  coefficients  are  given  in  (3.2.2).  The 
validity  of  this  formula  has  not  been  confirmed  by  experiment;  predictions  in  the  acoustic 
domain  are  typically  20  dB  smaller  than  those  given  by  Sevik’s  formula  (3.2.6)  when  u8/U  « 
20  -  30. 


The  value  of  e  («  0.2)  determines  the  magnitude  of  the  spectral  peak  in  Figure  3.1.1b  at 
k  =  |ac0|  when  M  1,  although  the  existence  of  this  theoretical  maximum  has  not  been 
confirmed  experimentally.  The  value  of  e  is  governed  by  the  size  of  the  boundary  layer, 
the  attenuation  and  refraction  of  surface  propagating  sound  waves  by  the  boundary  layer, 
by  diffraction  when  the  wall  has  a  finite  radius  of  curvature,  and  by  scattering  at  surface 
discontinuities  and  edges. 

§3.3  Boundary  layer  generated  sound 

At  low  Mach  numbers  the  wavelength  of  sound  generated  by  the  boundary  layer  greatly 
exceeds  its  thickness  S,  and  there  exists  a  simple  relation  between  the  acoustic  domain  of 
P(k,  u>)  and  the  frequency  spectrum  of  the  sound.  When  M  is  not  small  this  relation  is  com¬ 
plicated  by  refraction,  scattering  and  by  sound  generation  by  turbulence  in  the  outer  regions 
of  the  boundary  layer  at  distances  from  the  wall  comparable  to  the  acoustic  wavelength. 

Consider  the  sound  radiated  from  a  large  region  of  the  boundary  layer  of  area  A,  which  may 
be  assumed  to  be  centered  on  the  origin.  At  low  Mach  numbers  the  acoustic  wavelengths 
are  very  much  larger  than  the  boundary  layer  thickness  ( k08  <C  1),  and  the  pressures  at 
these  wavenumbers  are  constant  across  the  boundary  layer;  the  acoustic  pressure  frequency 
spectrum  $(x,  u)  of  the  boundary  layer  noise  is  then  given  by 

<E>(x,o>)  =  2^/C°  c°S-  -- P (ac0x/|x| ,  a;) ,  |x| 00,  where  (p2(x,t))  =  f  $(x,w)dw. 

|x|  Jo 

(3.3. 
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where  0  is  the  angle  between  the  radiation  direction  x  and  the  normal  to  the  wall  (the 
x2-axis). 

There  are  no  universally  accepted  measurements  of  the  wall  pressure  spectrum  in  the 
acoustic  domain  for  use  in  (3.3.1).  Sevik’s  wavenumber-white  approximation  (3.2.6)  yields 

$(x,o>)  ^  25iAM^v*/U)  cos2  0  3  g. 

plv%8  |x|2(o;5/{7)5 

according  to  which  the  peak  radiation  direction  is  normal  to  the  wall.  A  more  complicated 
dependence  on  radiation  direction  is  predicted  by  the  Chase  formula  (3.2.8).  Unless  u6/U  is 
very  large,  predictions  made  with  Sevik’s  formula  are  larger  (by  about  20  dB  at  u5/U  «  20 
for  M  =  0.01).  Thus,  equation  (3.3.2)  should  probably  be  regarded  as  an  upper  bound  for 
the  quadrupole  noise. 

§3.4  Rough  walls 

The  wall  pressure  wavenumber-frequency  spectrum  for  flow  over  a  rough  wall  whose  char¬ 
acteristic  roughness  height  R  is  much  smaller  than  5  can  be  approximated  by  a  smooth  wall 
empirical  formula  provided  the  friction  velocity  is  increased  to  compensate  for  the  increased 
surface  drag  and  turbulence  production  [30].  This  procedure  is  satisfactory  in  the  convec¬ 
tive  domain,  but  does  not  account  for  changes  at  low  wavenumbers  which  are  produced  by 
scattering  of  convective  pressures  by  roughness  elements.  In  particular,  it  seems  that  when 
uR/v,  >  5,  the  aerodynamic  sound  generated  by  the  boundary  layer  is  dominated  by  this 
scattering. 

To  investigate  the  influence  of  surface  roughness,  let  us  assume  that  the  motion  is  at 
low  Mach  number  and  of  sufficiently  high  Reynolds  number  that  viscous  stresses  can  be 
disregarded;  this  will  be  the  case  if  the  roughness  elements  protrude  beyond  the  viscous 
sublayer,  i.e.  for  Rv+/v  >  10. 

The  roughness  elements  behave  like  point  dipoles,  the  dipole  strength  being  the  unsteady 
drag  on  an  element.  If  they  are  modeled  as  rigid  hemispheres  distributed  randomly  over 
the  wall  (Figure  3.2a)  the  following  empirical  spectrum  is  consistent  with  measurements  by 
Hersh  [31]  of  the  noise  generated  by  low  Mach  number  flow  of  air  from  a  sand  roughened 
pipe  (Figure  3.2b),  provided  R  is  identified  with  the  mean  roughness  height: 


$ft(x,a;)  Acos2  9  Rvl  (uR/v*)3 
P20vl5  ~  [1  +  /?M/u,)2]5-5  ’ 
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*2 


Figure  3.2.  (a)  Wall  roughness  modelled  by  rigid  hemispheres, 
(b)  Roughness  noise  measured  by  Hersh  [31]. 
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where  r0  and  (3  are  empirical  constants,  and  9  is  the  angle  between  the  radiation  direction 
x  and  the  mean  flow.  The  best  fit  to  the  shape  of  the  experimental  spectrum  is  obtained 
for  (3  =  0.0025,  as  illustrated  by  the  solid  curve  in  Figure  3.2b,  whose  absolute  position 
has  been  adjusted  to  overlie  the  data,  since  no  directivity  information  is  available  from 
the  Hersh  experiments.  This  value  of  (3  and  the  exponent  5.5  of  the  denominator  in  (3.4.1) 
should  be  regarded  as  tentative.  The  “roughness”  parameter  t0  depends  on  the  spacing  of  the 
roughness  elements,  determined  by  the  fraction  a  of  the  plane  wall  covered  by  hemispheres. 
When  a  <  1,  it  appears  that  r0  «  {v*/Uc)2(a/ 7r). 

§3.5  Homogeneous  elastic  walls 

For  a  smooth,  homogenous  and  compliant  wall  the  boundary  layer  noise  at  very  small 
Mach  numbers  is  given  by  the  following  modification  of  equation  (3.3.1) 

4>(x,W)  =  ^|^^|1  +  R|2P(-c<,x/|x|.W),  |x[  -4  oo,  (3.5.1) 

z|x| 

where  R  is  the  “reflection  coefficient”  for  a  plane  acoustic  wave  of  frequency  u  specularly 
reflected  at  the  wall  into  the  radiation  direction  x. 

For  example,  when  it  is  permissible  to  model  the  wall  as  a  vacuum  backed ,  thin  elastic 
plate  of  bending  stiffness  B  and  mass  m  per  unit  area, 

_  (a;/a>c)  cos  9{(a;/a;c)2  sin4  ©  -  1}  +  ie  .  g  2x 

(u/uc)  cos  0{(a>/a>c)2  sin4  ©  —  1}  —  it 

where  u)c  =  <?0\Jrnj B  is  the  coincidence  frequency  (above  which  the  phase  speed  of  bending 
waves  on  the  plate  in  vacuo  exceeds  the  speed  of  sound  c0) ,  and  e  is  a  fluid  loading  parameter 
given  by 

_  PoC 1  (  1-2(7 

6  psc0  \12(1  -  or)2 

where  cr,  ps,  ci,  E  are  respectively  the  Poisson’s  ratio,  mass  density,  P-wave  speed,  and 
Young’s  modulus  of  the  material  of  the  plate  [32].  For  a  steel  plate  in  water  e  «  0.135;  for 
aluminum  in  water  e  ~  0.391.  For  an  aluminum  or  steel  plate  of  thickness  h  in  water  the 
coincidence  frequency  ojc  «  0.95 c0/h. 


E 


I2psc20{l  -  a2) t 


(3.5.3) 
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§3.6  Sound  generated  by  a  surface  irregularity 

Consider  low  Mach  number  turbulent  flow  past  a  surface  irregularity,  such  as  the  hemi¬ 
spherical  boss  in  Figure  3.3a  attached  to  a  plane  surface  x2  =  0.  The  acoustic  radiation 
can  be  calculated  from  (2.6.4)  when  the  vorticity  and  velocity  distributions  are  known.  The 
compact  Green’s  function  depends  on  the  wall  properties: 

Rigid  wall 

G(x'y’*-r)=4^^i(t-r-|x-Y|/c”)+spr^i(t-T— |x-yM;  (3.6.1) 


Elastic  wall 


G(x,y,t-  t)  =  |X- Y|/c„) 

+  8  /_“  R(o>) exp  (  -  iu,(t  -  r  -  |X  -  Y|/c„))  da,.  (3.6.2) 

In  these  formulae 

Y  =  (71,3/2,  n),  Y  =  (71,-3/2,73), 

where,  for  j  =  1  and  3,  Yj  is  the  velocity  potential  of  irrotational  flow  at  unit  speed  past  the 
surface  irregularity  in  directions  parallel  to  the  wall.  R(w)  is  the  wall  reflection  coefficient, 
defined  as  in  §3.5.  For  a  rigid  wall  R  =  1,  and  (3.6.2)  reduces  to  the  hard  wall  case  (3.6.1). 
For  a  vacuum-backed  thin  plate,  R  is  given  by  (3.5.2). 


The  radiation  is  of  dipole  type,  and  can  be  expressed  in  terms  of  the  unsteady  drag  force 
Di(t )  exerted  on  the  irregularity  (c.f.  (2.6.7)).  If  the  coordinate  origin  is  taken  within  the 
irregularity  we  find,  as  |x|  — >  oo: 


p(x,t) 


27TC, 


for  a  rigid  wall; 


|x|/c„ 


■Xj  \dDj(t)~\ 
:0|x|2  [  dt  Jt_ 


-Xi 


-|x|/Co 


=  - f°°  Di( u))R{u)e~iwt  du  for  an  elastic  wall,  (3.6.3) 

47TC0|x|2  [  dt  dtJ- oo  lt-\x\/c0 

where  Di(u)  =  (1/27t)  Di(t)eMt  dt.  In  the  elastic  case  this  formula  assumes  the  mass  of 
the  irregularity  is  small  compared  to  the  mass  of  the  area  of  the  plate  that  it  covers  (so  that 
it  has  a  negligible  influence  on  flexural  motions  of  the  plate) 
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Figure  3.3.  (a)  The  strength  of  the  dipole  radiation  from  a  small  surface 
irregularity  is  determined  by  the  drag  force  D{. 

(b)  Aerodynamically  smooth  plate  with  mass  load. 
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§3.7  Pinned  elastic  plate 

Consider  a  thin,  homogeneous,  vacuum-backed  elastic  plate  of  bending  stiffness  B  and 
mass  density  m  per  unit  area  loaded  with  a  point  mass  m0  at  x  =  0  (Figure  3.3b),  such 
that  the  surface  of  the  plate  exposed  to  flow  is  aerodynamically  smooth.  Green’s  function  is 

given  by  [33] 

m0  cos  0  rr  “3{  1  -  R(^)}ei(k'y+7(fc)y2"w[tl)  cPkdw 
G(x, y,  t  -  t)  =  G.(*.  y,  t  -  T)  -  2{2lr).  |X|  7(Jb)Z(u)D(t, u>) 

(3.7.1) 

where  G0(x,y,t  -  r)  is  defined  by  (3.6.2)  with  Y  =  Y  =  y  (for  a  smooth  wall).  In  the 
integral,  [t]  =  t-r  -  |x|/c0  denotes  retarded  time,  and 

7(&)  =  sgn {u)\f<?0  -  /c2,  +iy/k^~K 2  according  as  fc  >  |k0|; 


2?(fc,w) 


Bk4  —  mu2  - 

1  _  m°K°  [ 
2nm  Jo 


ipoQJ2  m 

7 (k)  ’ 

>  XdX 

[(1  -  iV)  A4  - 1  -  (e/pyv^^] 


K0  = 


the  vacuum  bending  wavenumber, 


(3.7.2) 


and  77  is  a  loss  factor  for  the  plate  (<  0  according  as  u  <  0).  When  w  >  0  the  square  root 
in  the  integral  for  2(u>)  is  negative  imaginary  when  A2  <  the  value  of  Z(ui)  for  u  <  0  is 
given  by 

Z(u>)  = 

where  the  asterisk  denotes  complex  conjugate. 


In  the  limit  m0  ->  00  the  plate  is  “pinned”  at  x  =  0,  and  cannot  execute  flexural  motions 
(although  bending  still  occurs  there).  This  might  be  taken  to  model  the  effect  of  an  isolated 
rivet;  if  the  “rivet  head”  is  not  flush  with  the  surface  exposed  to  the  flow  equation  (3.7.1) 
is  modified  by  taking  G0  to  be  defined  exactly  as  in  (3.6.2),  with  Y,  Y  determined  by  the 
shape  of  the  rivet  head.  In  this  case  there  are  two  components  of  the  sound  generated  at 
the  rivet,  corresponding  respectively  to  the  drag  dipole  (3.6.3)  and  the  pinning  contribution 
given  by  the  integral  term  in  (3.7.1). 

When  the  radiation  is  expressed  in  terms  of  the  blocked  pressure  ps  we  make  use  of  the 
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reciprocal  theorem  formula  (2.7.3),  in  which 


8Cr  _  (dG\ 

8t  ~{dy2)y2= o' 


The  contribution  to  the  far  field  sound  from  mass  load  m0  (i.e.,  from  the  final  term  in  (3.7.1)) 
can  then  be  expressed  in  the  form 

— zm0  cos  9  r  r 00  oj3{1  -  R(a;)}pJi(k,a;)e~t'a,(t~|x|/Co)  j2, 

87r2c0|x|  JJ-oo  Z(u)T>(k,u) 


p(x,  t) 


d2kdu, 


oo. 


(3.7.3) 


This  integral  will  have  significant  contributions  from  two  distinct  wavenumber  regions:  (i) 
from  the  convective  range  of  wavenumbers  (k\  ~  u/Uc,  k$  ~  0)  that  defines  the  dominant 
wall  pressure  fluctuations,  and  (ii)  from  the  vicinity  of  the  region  where  V(k,u)  =  0,  which 
corresponds  to  the  sound  radiated  when  bending  waves  generated  by  subconvective  wall 
pressures  over  the  whole  plate  are  scattered  by  the  mass  load.  The  actual  contribution  in 
case  (ii)  will  be  determined  by  the  size  of  the  boundary  layer  flow  and  the  magnitude  of  the 
loss  factor  77. 


When  m0  ->  00  (plate  pinned  at  x  =  0),  the  acoustic  pressure  is  independent  of  m0,  and 
can  be  written 

im cos 0  h  r  r°°  oj|i<;|{1  —  R(a;)}ps(k,a;)e~t^t~lx^Co^  ^ 

47t(wc/i/c0)  |x|  J  /-oo  l(u>)V(k,L)) 


p(x,  t) 


crkdu;,  |x|  — >  00, 

(3.7.4) 


where 


XdX 


[(1  -  irj) A4  -  1  -  (e/p)/^X2  -  n2 


§3.8  Forward  and  backward  facing  steps 

Wall  surface  irregularities  in  the  form  of  forward  or  backward  facing  steps  produce  drag 
fluctuations  that  are  equivalent  to  localized  aeroacoustic  dipole  sources.  To  fix  ideas,  con¬ 
sider  the  forward  facing  step  of  height  h  illustrated  schematically  in  Figure  3.4.  Take  the 
coordinate  origin  at  the  foot  O  of  the  step,  with  the  mean  flow  locally  in  the  aq-direction. 
The  sound  radiated  by  low  Mach  number  turbulent  flow  over  the  step  is  calculated  using  the 
Green’s  function  (3.6.1)  when  the  surface  is  rigid.  When  the  step  is  formed  by  two  overlap¬ 
ping,  identical  plates,  the  step  noise  is  determined  by  (3.6.2),  where  R(w)  is  the  reflection 
coefficient  for  a  plane  wave  incident  on  an  infinite  homogeneous  plate.  Because  the  wall  is 
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uniform  in  the  z3-direction  we  have  Y3  =  y3;  the  function  Yx  is  determined  by  conformal 
mapping  of  the  complex  region  z  —  yx  +  iy2  “above”  the  step  onto  the  upper  half  Im  (  >  0 
of  the  (-plane  by  means  of  the  formula 

f  *=  i  =F  In  (<  +  ±  1),  ta<>0,  (3.8.1) 

where  the  upper/lower  signs  are  to  be  taken  respectively  for  forward/backward  facing  steps. 
The  foot  of  the  step  0  maps  into  the  point  (  =  TY  and  the  top  S  maps  into  (  =  ±1.  Then 

Y,(y)  =  R*  (f)  •  (3-8.2) 


When  the  step  height  h  is  small  compared  to  the  dominant  turbulence  length  scale  we  can 
take 


h  roo  g*(fc!/l+7(%2) 

^i(y)  -  vi  ow  / 


ih 


(3.8.3) 


The  Green’s  function  can  be  used  if  the  vorticity  and  velocity  distributions  are  known.  If, 
however,  the  unsteady  drag  Dx(x3,  t),  per  unit  span  (^-direction)  can  be  found  by  alternative 
means,  (3.6.3)  gives  the  radiation  in  the  form, 


p(x,t)  ~ 


—  cos#  d 
47rc0|x|  dt 

—  cos#  d 
47tc0|x|  dt 


f  dy3  f  Dx(y3,u){l  +  R(u>)}e  du), 

J  span  J  — oo 

[  {£>,(»,[*])+  r  Di(</3,w)R(w)e-<"Mrfw}  dy3, 

J  span  L  J— oo  ) 


|x|  — »■  oo, 
(3.8.4) 


where  6  is  the  angle  between  the  radiation  direction  x  and  the  flow  direction  (xi-axis),  and 
[f]  =  t  -  \Jx\  +  x%  +  (x3  -  yz)2/c0.  R  =  1  for  a  rigid  wall,  and  then 


P(X' ()  “  2^\ItLn  Dl (V 3>  W)  dV3'  W 


27rc0|x|  dt  J  span 


(3.8.5) 
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§3.9  Rib-stiffened  elastic  wall 

Suppose  the  wall  can  be  modeled  as  an  aerodynamically  smooth,  vacuum-backed,  thin 
elastic  plate  with  a  narrow  rib-stiffener  along  the  x3-axis  (Figure  3.5).  The  sound  radiated 
from  the  rib  can  be  calculated  from  the  reciprocal  theorem  formula  (2.7.3)  in  terms  of  the 
blocked  wall  pressure  ps. 

When  the  dominant  wall  pressure  fluctuations  have  length  scales  larger  than  the  width  of 
the  rib,  the  far  field  acoustic  pressure  at  x  can  be  expressed  in  the  general  form  [34] 

7r m  d  rr°°  u;2{l  -  R(o;)}aj(ni,^)(-fe/K0)Jp,(fc,/c0n3,a;)e~^(t~|x|/c°)(ifcda> 
P(X,f)  “  K'Vy#  +  I$vi,w) 

(3.9.1) 

In  this  formula  ps(k,K0n3,u )  is  the  Fourier  transform  of  the  blocked  wall  pressure,  R(w) 
is  the  reflection  coefficient  (3.5.2),  and  n  =  (ni,n2,n3)  =  x/|x|  is  the  unit  vector  in  the 
radiation  direction. 


The  coefficients  aj(n uu)  are  dimensionless  and  depend  on  the  way  in  which  the  plate  is 
attached  to  the  rib.  When  the  plate  is  clamped  to  the  rib  (so  that  (  =  0  and  dQ/dxi  =  0  at 
arj  =  0,  C  =  plate  flexural  displacement)  we  have  [34] 


^2  +  ^ini/i  -[^i  +  ^onip] 

a°~  Wi  ’  ai  “  -  Wi 


a2  =  a3  =  0, 


(3.9.2) 


where 

*;  = 

For  w  >  Owe  take  \/A2  —  JP  >  0  when  |A|  >  —  n3,  and  \M2  —  p2  —  —iy/p?  —  A2  for 

|A|  <  pyfl  —  n3.  For  <  0  the  ctj  are  determined  from  the  formula 

ajfau-u)  =  (-lya^nuu). 

To  calculate  the  acoustic  pressure  frequency  spectrum  (defined  as  in  the  second  of  equa¬ 
tions  (3.3.1))  for  a  rib  of  wetted  span  L  we  use  the  result 


/. 


\jd\ 


A2  =  A2  +  p2n\. 


-oo  (1  -  ir])  A4  -  1  -  (e/p)/y/7?  -  p2 


(3.9.3) 


{ps{ki,k3,aj)p*(k[,k3,u}'))  «  —6(u>  -  u')6(ki  -  k[)P(ki,k3,u),  L~>8.  (3.9.4) 

Z7T 
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Then  (3.9.1)  supplies 


7tLm2u5\l  —  R(cu)l2 


P  {k,K0n3,u) 

EJ-aj(ni,a;)(- 

k/K.)’ 

2 

dk 

V{ 

fsJk2  +  Klnl,w) 

2 

|x|  — >  oo,  u)  >  0.  (3.9.5) 


The  remarks  following  equation  (3.7.3)  concerning  the  dominant  contributions  to  the  in¬ 
tegral  in  (3.7.3)  also  apply  to  equations  (3.9.1)  and  (3.9.5),  namely,  that  the  integrals  will 
tend  to  be  dominated  by  contributions  from  (i)  from  the  convective  region  ki  ~  w/UC)  and 
(ii)  resonant  wall  modes  where  f,u/)  <  1,  the  latter  being  controlled  by  the 

boundary  layer  size  and  the  loss  factor  rj. 
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§3.10  Flow-excited  panel  radiation 

The  method  of  §3.9  is  easily  extended  to  the  case  illustrated  in  Figure  3.6,  in  which  an 
elastic  plate  is  divided  into  equal  strip-like  panels  transverse  to  the  mean  flow  direction  [34], 
the  ribs  being  located  at  xx  =  Nd,  -oo  <  N  <  oo.  The  acoustic  pressure  may  now  be 
expressed  in  the  form 

~  2^1  yv  y  f°°  a;2{l  -  R(a>)}(-l)Jaj(ni,q;)  //c0nx  +  2 Nn[d\3 

P  c0d|x|  dt  N=_ooj=0  J- oo  K0V{^J (K0ni  +  2A^7r/d)2  +  w)  \  K°  ) 

X  p.(<e»i*i  +  |x| ->  oo.  (3.10.1) 


For  clamped  panel  edges  the  coefficients  olj  are  given  by  equations  (3.9.2),  but  with  tyj 


defined  by 


2it  _ A# _ 

K°d  n=- oo  (1  -  iri)h.%  -  1  -  (e/n)/^Ah  -  /i2 


(3.10.2) 


where 

Ajv  =  2Ntt/ K0d  —  fin i,  =  A^  -I-  fj^n^. 


The  formulae  (3.10.1)  and  (3.10.2)  are  exact  within  the  approximation  of  linear,  thin  plate 
theory,  and  with  the  neglect  of  mean  flow  convection  of  sound.  However,  correspondingly 
simple  formulae  are  not  available  for  panels  of  finite  dimensions,  such  as  the  rectangular 
panels  relevant  in  underwater  acoustics.  In  this  case  the  reciprocal  formula  (2.7.3)  gives  a 
formally  exact  procedure  for  calculating  the  radiation  in  terms  of  the  wall  pressure  fluctua¬ 
tions,  and  the  flexural  displacement  (R(x,  y,  u)  of  an  arbitrary  point  y  of  the  panel  produced 
by  a  time  harmonic  volume  source  of  unit  strength  placed  at  the  far  field  point  x.  This  re¬ 
duces  the  problem  to  the  solution  of  a  well  defined  diffraction  problem,  which  must  be  solved 
numerically. 
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Figure  3.6.  Periodically  rib-stiffened  wall. 
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4.  SOUND  GENERATED  AT  THE  EDGE  OF  A  LARGE  SURFACE 

In  underwater  applications  free  stream  turbulence  inflow  on  a  leading  edge  is  not  usually 
an  important  source  of  sound,  except  possibly  when  the  impinging  vorticity  is  shed  from  a 
neighboring  upstream  solid  boundary  or  control  surface.  The  contribution  from  trailing  edges 
also  tends  to  be  small  because  the  upwash  velocity  generated  by  inflow  turbulence  is  largely 
canceled  by  shedding  from  the  edge.  Nevertheless,  experiments  at  high  Reynolds  numbers 
[35  -  37]  indicate  that  a  trailing  edge  can  be  an  important  source  of  high  frequency  “self¬ 
noise”  ,  attributable  to  relatively  small  scale  turbulence,  generated  by  the  natural  instability 
of  the  boundary  layers  on  the  airfoil,  interacting  with  the  trailing  edge.  At  lower  Reynolds 
numbers  (smaller  than  about  2  x  105  when  based  on  airfoil  chord),  quasi-periodic  shedding 
of  vorticity  can  occur,  especially  when  the  trailing  edge  is  blunt,  causing  the  airfoil  to  “sing” 
at  a  frequency  uj  «  U/h,  where  h  is  the  airfoil  thickness  near  the  edge  and  U  the  mean 
stream  velocity  [30]. 

§4.1  Edge  noise  Green’s  functions 

At  very  small  Mach  numbers  the  wavelength  of  the  sound  produced  by  the  interaction  of 
unsteady  flow  with  the  edge  is  always  very  much  larger  than  the  “hydrodynamic  wavelength”, 
which  is  effectively  the  distance  between  the  edge  and  those  components  of  the  turbulence 
that  interact  effectively  with  the  edge.  This  makes  it  possible  to  express  the  edge  noise  in 
terms  of  a  compact  Green’s  function,  defined  such  that  the  source  point  y  is  very  much  closer 
to  the  edge  than  a  characteristic  acoustic  wavelength.  We  give  expressions  for  the  frequency 
domain  compact  Green’s  function  G(x,  y, ;  uj)  for  the  semi-infinite  configurations  illustrated 
in  Figure  4.1: 


i.  Rigid  half-plane  [11]  (Figure  4.1a): 

The  plane  occupies  x\  <  0,  Xz  =  0 

G(x,  y;  uj)  =  G0(x,  y;  uj)  +  Gx  (x,  y;  uj)  + . . . , 

where,  for  |x  -  y3i3|  ->  oo  and  k0\/2/i +2/1  <  1, 


(4.1.1) 


Go(x,y;w)  = 


-1 


47t|x  -  y3i3 


.^clx-yaM  Gi  (x,  y ;  uj) 


-1  V^^*(x)y?*(y)  fejx-Mhi 

'Ksfisiri  |x-y3i3|3/2 

(4.1.2) 
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semi-infinite  rigid  plate 


•  x 


compact  thickness  with  rounded  edge 

•  x 


Figure  4.1.  Edge  scattering  configurations. 


Section  4.1 


52 


Large  Surface  with  an  Edge 


Report  No.  AM  98-001 


Boston  University,  College  of  Engineering 


i3  is  a  unit  vector  parallel  to  the  rr3-axis  (the  edge),  and 

ip*(x )  =  y/rsin(6/2)  (4.1.3) 

is  a  velocity  potential  of  incompressible  flow  around  the  edge  (in  the  anticlockwise  direction) 
expressed  in  terms  of  polar  coordinates  (xi,x2)  =  r(cos0,sin0).  The  component  G0  of  G  is 
the  radiation  from  a  point  source  at  y  when  scattering  is  neglected.  Gi  is  the  first  correction 
due  to  the  presence  of  the  half-plane,  and  may  also  be  expressed  in  the  integral  form 

sgn(x2)sgn(y2)  e,'lfcyi+fclI1+fc3(x3"y3)+^|y2|+7|:C2|}dA:d/c1dfc3 
Gx(x,y;o;)  =  +  ^  +  ¥l(k  Th+  i0) 


—m*  (y)sgn(a:2)ei?  ff°°  ei^1I1+7|l2|+/C3(x3-J'3)}  - r  , 

-  —  V  -- — //  - /_  ,  , - dkidhz,  K’oyVi  +  y2  <  1, 

47T2  JJ-oo  V^o  +  rui 

(4.1.4) 


where  _ 

7  =  \JKl  -  fe?  ~  7  =  \Jk2  —  k2  —  k%,  R0  =  y/Kl-kl 

are  either  real  with  sign  sgn(u;)  or  positive  imaginary. 

The  component  G\  of  the  compact  Green’s  function  reduces  to  a  particularly  simple  form 
for  two-dimensional  problems  involving  sources  with  no  dependence  on  the  spanwise  coordi¬ 
nate  x3.  The  dependence  of  Gi  on  x3,  y3  is  eliminated  by  integrating  over  -oo  <  y3  <  oo, 
using  the  method  of  stationary  phase  for  K0yjx \  +  x\  — >  oo.  When  the  result  is  multiplied 
by  (— l/27r)e-icJ^-T^  and  integrated  over  all  frequencies,  we  find 

Gi(x,y,t-r)  »  -  ^-6(t  —  r  —  |x|/c0),  |x|  -»  oo,  (4.1.5) 

where  x  =  {x\,x2),  y  =  (yi,  y2)  in  two-dimensions.  The  two-dimensional  counterpart  of  Go 
is  independent  of  y  and  is  therefore  of  use  only  in  applications  involving  monopole  sources, 
which  are  uninfluenced  by  the  edge  to  first  order. 

ii.  Rigid  half-plane  with  attachment  (Figure  4.1b): 

The  representations  (4.1.1)  -  (4.1.4)  remain  valid  when  the  edge  of  the  half  plane  supports 
a  rigid,  “attachment”  whose  cross-section  is  acoustically  compact  and  length  scale  of  varia¬ 
tion  in  the  spanwise  (x3)  direction  is  compact,  provided  the  potential  <p*( y)  =  </?*(yi,  y2)  is 
replaced  by  <f>*(y)  =  <h*(yi,  y2,  y3)  which  describes  potential  flow  around  the  irregular  edge, 
such  that 

S/2, 2/3)  -*¥>*(i/i,ife),  as  \jyl  +  yl  ->  00. 
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iii.  Rigid  half-plane  of  compact  thickness  (Figure  4.1c): 

Formulae  (4.1.1)  -  (4.1.4)  can  be  used  provided  the  wall  thickness  h  is  acoustically  compact 
and  <p*(y)  is  re-defined  as  in  Case  ii.  For  example,  for  the  square-edged,  thick  half-plane 
of  Figure  4.1c  (where  the  negative  aq-axis  is  assumed  to  lie  in  the  upper  surface,  with  the 
origin  at  the  edge)  we  find  by  conformal  transformation 

f(y)  = 


|  =  /(C)  £  - ifCi/C2  -  1  -  >n«  +  %/C2  -  1)}  -  *.  (4-1.6) 

where  the  transformation  in  the  second  line  maps  the  fluid  region  in  Figure  4.1c  onto  the 
upper  half  (Im(  >  0)  of  the  (-plane,  where  VC2  -  1  is  taken  to  have  positive  imaginary  part. 

iv.  Rigid  wedge  (Figure  4. Id): 

The  dominant  term  in  the  compact  Green’s  function  for  a  rigid  wedge  of  interior  angle  x 
can  be  expressed  in  the  following  manner,  analogous  to  Gi(x,y;cu)  of  (4.1.2): 


-1 


«£p*(x)y*(y) 


gfrIVQ  1  2/3l3| 


2/3*31  “>  00,  V  = 


G(x,y;w) «  |x  -  „3i3|>+- 

where  i3  is  a  unit  vector  in  the  rc3-direction  (parallel  to  the  edge),  and 


7T 


2x-x 


,  (4.1.7) 


V?*(x)  =  Re^e-^),  z  =  xx  +  ix2,  (4.1.8) 

is  a  velocity  potential  of  incompressible  flow  around  the  edge  (in  the  anticlockwise  direc¬ 
tion).  When  the  wedge  is  rounded,  or  has  an  acoustically  compact  attachment,  of  the 
type  considered  for  the  half  plane,  the  potential  <p*(y)  =  ¥>*(2/1, 2/2)  must  be  replaced  by 
$*(y)  =  $*(yi,  y2, 2/3)  which  describes  potential  flow  around  the  irregular  edge,  with  the 
asymptotic  property 

$*(l/i,  !fe,  2/3)  ->  Re^e-^) ,  z'  =  yi  +  iy2,  as  \Jy\  +  yl  -4  00.  (4.1.9) 


For  the  half-plane,  the  component  Go  of  G  represents  the  radiation  from  a  point  source  at 
y  when  scattering  is  neglected;  Gi  is  the  first  correction  due  to  presence  of  the  half-plane. 
Thus,  for  the  Reynolds  stress  fluctuation  (p0u;Wj)(x,  w)  near  the  edge,  the  edge  generated 
acoustic  pressure  is  given  by 


p{x,u) 


J  (PoViVj)(y,u) 


dV(y)  y*(x) 
dyidyj  |x-y3i3|3/2 


exp  (z«0|x  -  y3i3|)  d3 y,  |x-y3i3|  ->  00. 
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In  the  time  domain  p{x,t )  =  f^p(x,u)e~iu,tdu;,  but  it  is  not  necessary  to  evaluate  the 
integral  to  estimate  the  acoustic  efficiency.  For  an  eddy  of  dimension  i  and  frequency  v/£, 
distance  4  from  the  edge,  we  find 

£  / £\ f  _  ,  X 

p(x,  t)  ~  1^-1  p0v2VM sin5  ip sin(-0),  M  ~  v/c0  <1,  |x|  oo, 

where  ip  is  the  angle  between  the  observer  direction  and  the  edge  of  the  plate.  The  acoustic 
pressure  is  a  maximum  in  the  upstream  directions  0  ~  ±7r, 

Thus,  the  acoustic  efficiency  ~  ( £/£0)3M 2,  which  is  larger  than  the  efficiency  of  turbulence 
near  a  compact  solid  by  a  factor  1/M.  For  the  case  of  turbulence  in  a  boundary  layer  flow 
over  the  edge,  for  example,  the  edge-generated  sound  power  is  the  same  as  the  direct  acoustic 
power  produced  by  turbulence  quadrupoles  (of  efficiency  ~  M5)  in  a  region  of  the  boundary 
layer  extending  a  distance  4/M3  »  £q  upstream  of  the  edge.  In  water  this  distance  can  be 
so  large  that  all  of  the  noise  can  be  considered  to  be  generated  at  the  edge.  The  frequency 
of  the  edge  noise  will  also  tend  to  be  higher  by  a  factor  ~  U/v,  where  U  and  v  are  the  mean 
and  fluctuating  velocity  components,  because  the  frequency  of  the  quadrupole  sound  ~  v/S, 
where  6  is  of  the  order  of  the  boundary  layer  thickness,  whereas  for  edge  noise  the  frequency 
~  U/5  is  set  by  the  rate  at  which  the  turbulence  is  swept  past  the  edge,  where  d2ip* /dyidyj 
is  large. 

§4.2  General  representation  of  trailing  edge  noise  [38  -  41] 

Consider  turbulent  flow  past  the  trailing  edge  of  the  upper  wall  (x2  =  0)  of  a  non-compact 
surface,  which  may  be  regarded  as  any  one  of  the  surfaces  in  Figure  4.1.  The  main  stream 
outside  the  boundary  layer  has  low  subsonic  speed  U  in  the  rrq-direction  (as  indicated  in 
Figure  4.2  for  a  thin  half-plane).  The  details  of  the  boundary  layer  turbulence  are  unknown, 
but  it  will  be  assumed  that  during  the  interaction  with  the  edge  the  statistical  properties  of 
the  impinging  turbulence  are  the  same  as  those  several  boundary  layer  thicknesses  upstream. 
This  permits  the  edge  noise  to  be  expressed  in  terms  of  the  boundary  layer  wall  pressure  in 
the  upstream  region,  which  may  regarded  as  uninfluenced  by  the  proximity  of  the  edge. 

The  calculation  of  the  noise  can  be  formulated  as  a  scattering  problem,  in  which  the 
pressure  p\,  say,  that  would  be  produced  by  the  same  turbulent  flow  if  the  surface  were 
removed,  is  scattered  by  the  edge.  The  scattered  pressure  p'  includes  both  acoustic  and 
hydrodynamic  components,  the  latter  accounting  for  the  modification  of  the  near  field 
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Figure  4.2.  Coordinates  for  trailing  edge  noise. 
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pressure  by  the  surface.  The  condition  that  the  normal  velocity  vanishes  on  the  surface  will 
be  taken  in  the  high  Reynolds  number  approximation 

dpi/dxn  +  dp'/dxn  =  0,  on  S,  (4.2.1) 

where  S  denotes  the  entire  surface  (including  the  “upper”,  “lower”  and  “edge  region”). 
We  assume  that  in  turbulence-free  regions  pressure  fluctuations  of  frequency  u  satisfy  the 
Helmholtz  equation 

(V2  +  k20)P  =  0.  (4.2.2) 

No  explicit  account  is  taken  of  the  boundary  layer  turbulence,  except  insofar  as  it  is  responsi¬ 
ble  for  the  pressure  pp  in  particular  p'(x,  u)  is  assumed  to  be  a  solution  of  (4.2.2)  everywhere, 
so  that  the  scattering  of  sound  by  the  shear  flow  is  ignored.  This  approximation  is  not  valid 
at  very  high  frequencies  when  the  acoustic  wavelength  is  comparable  to  the  thickness  of  the 
boundary  layer  [42  -  44]. 

Because  of  its  definition,  the  pressure  pi(x,  u)  must  be  an  outgoing  solution  of  Helmholtz’s 
equation  in  x2  <  0,  in  the  region  “below”  the  boundary  layer  sources,  and  is  equal  to  half 
the  boundary  layer  blocked  pressure  ps  on  X2  —  +0  (c.f.  equation  (3.2.7)).  Hence, 

PiM  =  k,w)ei(kx"y(fc)X2)dA:1dfc3,  *2  <  0,  (4.2.3) 

where  7 (k)  is  defined  as  in  (4.1.4).  The  problem  of  determining  p'  now  reduces  to  a  consid¬ 
eration  of  the  scattering  of  each  Fourier  component  lp5(k,o;)e^k  x_7^^2^  of  p\  by  S. 

The  scattered  pressure  p'  is  an  outgoing  disturbance  both  above  and  below  S,  and  has  the 
following  general  representation  in  terms  of  the  Kirchhoff  formula 

p'(x,u>)  =  -  <f  y3,u>)C?(x,  y;  w)dS(y),  (4.2.4) 

Js  dyn 

provided  Green’s  function  is  chosen  to  have  vanishing  normal  derivative  on  S,  and  yn  is  normal 
to  S  and  directed  into  the  fluid.  At  low  Mach  numbers  the  distances  of  those  turbulence 
eddies  from  the  edge  that  are  responsible  for  the  edge  noise  are  always  very  much  smaller 
than  the  acoustic  wavelength,  so  that  G(x,y;u>)  may  be  approximated  by  the  generalized 
form  of  the  compact  Green’s  function  (4.1.2),  for  which  <p* (y)  is  replaced  by  $*(y)  (as 
described  in  §4.1ii,  iii). 
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When  the  result  of  this  calculation  is  integrated  over  all  frequencies,  the  pressure  in  the 
acoustic  far  field  is  found  to  be  given  by 


e-iu,«-|x| /co)dkldU),  |x|  oo. 

(4-2.5) 

where  9  =  axcsmfa/ \J x\  +  x$),  and  ip  is  the  angle  between  the  observer  direction  and  the 
edge  (the  z3-axis,  see  Figure  4.2  for  the  thin  half-plane).  The  peak  acoustic  pressures  are 
radiated  in  the  “forward”  directions  9  =  ±n.  The  influence  of  edge  geometry  is  governed  by 
the  complex  valued  function  1(h)  which  is  defined  by 

i(h)  =  c.c.  l^^e~iklZdz)  >  *i>o.  z  =  vi  +  m-, 


1(h)  =  c.c.  {i(-h)} ,  h<0 


(4.2.6) 


where  w(z)  is  the  analytic  function  whose  real  part  is  4>*(y),  “c.c.”  denotes  complex  con¬ 
jugate,  and  the  integration  is  taken  in  the  anticlockwise  sense  around  the  contour  C  formed 
by  the  surface  profile  of  S  in  the  y^-plane  .  If  a  spanwise  section  of  the  edge  of  length  L  is 
wetted  by  the  turbulent  flow,  the  relation  (3.9.4)  may  now  be  used  to  express  the  acoustic 
pressure  frequency  spectrum  in  the  form 


$(x,w) 


uL  sin2  (9/ 2)  sin  ip 
27r2c0|x|2 


f  \l(h)\2P(h,K>oCOsip,cj)  dh,  |x| ->  oo,  (4.2.7) 
J -00 


This  formula  is  applied  using  measured  values  of  the  wall  pressure  spectrum  P(k,u;).  In 
practice,  most  estimates  and  comparisons  with  experiment  are  based  on  a  measurement 
of  the  surface  pressure  frequency  spectrum  $pp(a;)  at  a  single  point,  which  is  then  used 
in  an  empirical  model  of  P(k,u;)  which  depends  linearly  on  <3>pp(u;).  The  simplest  is  the 
Corcos  model  (3.2.3),  which  represents  the  behavior  of  P(k,u;)  at  low  Mach  numbers  when 
oj6*/U  >  0.1  (5*  «  6/8)  in  the  immediate  vicinity  of  the  convective  ridge  (Figure  3.1),  where 
the  spectrum  has  a  large  maximum.  Because  of  the  dominant  contribution  from  this  region, 
the  integral  (4.2.7)  may  often  be  approximated  by  setting  h  =  u/Uc  in  1(h),  he- 


$(x,w) 


ojL  sin2  (9/2)  sin  ip 
27t2cJx|2 


/  (J  \  2  rOO 

l[jj  )  J  p (h,Ko  cos  ip,  v)dku 


X 


oo, 


(4.2.8) 
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§4.3  Influence  of  surface  geometry  on  edge  noise 

The  details  of  edge  geometry  become  important  at  high  frequencies,  when  the  turbulence 
length  scale  is  comparable  to  the  radius  of  curvature  of  the  edge.  The  following  examples 
illustrate  the  trends  to  be  expected  for  three  different  edge  types: 

Thin  half-plane  This  is  the  canonical  case  (a)  of  Figure  4.1,  for  which 

w(z)  =  T(fci)  =  —  e  4 )  ki  >  0.  (4.3.1) 


Square-edge  The  first  of  case  (c)  in  Figure  4.1;  compact  wall  thickness  h: 


/F  0.71  {1  +  0.2bkih  +  (kih)2} 

mkl)l  K  V*T(wjT - HW 


hh  >  0.1,  (4.3.2) 


where  £(z)  is  defined  in  (4.1.6) 

Rounded  edge  Probably  the  simplest  case  amenable  to  analytical  treatment  is  shown  in 
Figure  4.3b.  The  boundary  of  the  rounded  edge  is  mapped  onto  the  real  C-axis,  and  the  fluid 
is  mapped  into  Im  C  >  0  by 


z 

h 


/(C)  +  a/(  C/P) 

1  +  a 


and 


w{z ) 


h{l  +  a//32). 
tt(1  +  o) 


(4.3.3) 


where  /(C)  is  given  in  (4.1.6),  and  a  and  (3  are  numerical  coefficients.  For  the  case  shown  in 
Figure  4.3b,  a  =  0.0013,  (3  =  0.02,  the  mean  radius  of  curvature  of  the  edge  ~  l.lAh,  and 

\l(kl)\  «  ^  exp  -  ^-)  .  h  >  0.  (4.3.4) 

Half-plane  with  cylindrical  attachment  (Figure  4.3a).  The  semi-infinite  plane  (zi  < 
0,  X2  =  0)  is  smoothly  connected  at  its  edge  to  a  circular  cylinder  of  diameter  D  on  the 
underside  (x2  =  —0): 


w(z)  =  —  \/2xDC, 


ki  >  0, 


where  the  fluid  region  is  mapped  onto  the  upper  half  of  the  C-plane  by 

2  _  7T  C 

5”  l+Cln(l-l)  ' 


(4.3.5) 


(4.3.6) 
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*2 


0.1  1  10 


Figure  4.3.  Predicted  edge  noise  spectrum  (U/6*)$(x,u)/[a0(p0vt)2M{L5*/\-x.\2) sin  sin2  (0/2)] 
for  the  half-plane;  for  square  edge  of  thickness  h  =  <5;  for  a  half-plane  with  rounded 
attachment  of  diameter  D  =  and  for  a  rounded  edge  of  thickness  h  =  6. 
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When  the  Corcos  model  (3.2.3),  (3.2.4)  is  used,  the  acoustic  pressure  frequency  spectrum 
(4.2.8)  can  be  expressed  in  the  non-dimensional  form 


(C7&)*(x,w)  _ 

M)5  M 


'LL 


siin&sin2  (?)  (^)|l(|)|  j- 


aJjjjL/U) 

i + w/m 


3  ) 


a0  =  0.035, 
(4.3.7) 


where  M  =  U/c0  is  the  Mach  number  of  the  main  stream.  In  particular,  for  the  half-plane 


(C7/<?,)$(x,  u) 
{PoVlf 


w  M 


sin  ip  sin 


2 


aJjjL/U) 

[a%  +  (u5,/U)2]*  ‘ 


(4.3.8) 


The  general  expression  (4.3.7)  for  arbitrary  edge  geometry  reduces  to  this  result  when  the 
hydrodynamic  length  scale  UJoj  is  much  larger  than  the  characteristic  dimension  of  the  edge. 

The  predicted  forms  of  the  acoustic  pressure  spectra  are  compared  in  Figure  4.3  for  the 
above  cases  in  the  particular  case  in  which  h  =  D  =  S,  where  5  is  the  turbulent  boundary 
layer  thickness  for  flow  over  the  upper  surface. 

The  surface  pressure.  The  surface  pressure  on  the  upper  surface  of  a  thin  half-plane 
determined  by  the  scattering  theory  of  §4.2  can  be  cast  in  the  form 


p(Xl,±0,x3,u)  =  i  fr  p.( k,  w)  [l  ±  erffe-* -  fcfjv2  +  *,  )1  e‘k*  cPk,  x,  <  0, 
2"-“  1  J  (4.3.9) 

where  k  =  (ki,k3)  and  the  ±  sign  is  taken  on  the  upper/lower  surface  x2  =  ±0.  The 
argument  of  the  error  function  has  positive  real  part  for  all  real  values  of  k\,  so  that  the 
error  function  «  1  as  — >  -oo.  Thus,  p  -¥  0  on  the  lower  surface  (x2  =  -0)  far  upstream 

of  the  edge,  whereas  p->ps  =  2pi  on  the  surface  x2  =  +0  exposed  to  the  turbulent  stream. 
This  occurs  at  distances  upstream  of  the  edge  exceeding  the  characteristic  eddy  dimension, 
and  represents  a  pressure  doubling  by  “specular  reflection”  of  the  incident  pressure  p\  by  the 
plane,  rigid  surface.  Thus,  if  the  impinging  boundary  layer  turbulence  is  assumed  to  convect 
as  a  frozen  pattern  over  the  edge,  measurements  of  the  upstream  wall  pressure  ps  (or,  rather, 
of  its  wavenumber-frequency  spectrum)  can  be  used  to  predict  the  edge  noise  spectrum. 
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§4.4  Noise  from  flow  over  a  right-angled  wedge 

Consider  the  edge  noise  generated  by  turbulent  boundary  layer  flow  off  the  upper  surface 
x2  =  +0  of  the  rigid,  right  angled  wedge  bounded  by  the  planes  xi  <  0,  x2  =  0  and 
x\  =  0,  x2  <  0.  We  regard  the  turbulence  as  frozen  during  its  convection  past  the  edge, 
and  determine  the  sound  generated  for  the  two  wedges  illustrated  in  Figure  4.4:  (a)  with  a 
sharp,  and  (b)  with  a  rounded  edge. 

For  a  wedge  of  arbitrary  interior  angle  Xi  we  find  by  the  method  of  §4.2,  using  Green’s 
function  (4.1.7),  with  <p*(y)  replaced  by  <f>*(y)  for  a  rounded  edge  geometry,  that  the  acoustic 
pressure  is  given  by 


‘>  «  III  (*"  if’ ")  W  -  00, 

(4-4.1) 

where  v  =  7r/(27r  -  x),  and  1(h)  is  defined  as  in  (4.2.6)  with  the  path  of  integration  in  the 
z-plane  being  over  the  profile  of  the  wedge  in  the  yit/2-plane  in  the  anticlockwise  sense.  The 
corresponding  frequency  pressure  spectrum  for  a  wetted  span  of  length  L  is 


$(x,w) 


L  (  cos  v(t 
7r  \  r (v 

L  ( COS  v(6  —  7T 

7T  V  r(l/)|x| 


(0-7r)\  (u)  sin'll) 

)|x|  J  V  2c0  , 


2  v 
l 

/  \l(h)\2P(h,  K'o cos dh 

*  J—oo 


wsm 

~2ct 


p— ^  J  P(h,  Ko  COS  ip,  u)  dh  |x|  ->  oo. 

(4.4.2) 


For  a  sharp  cornered  wedge 


w  =  zv 


1(h) 


— 2v  sin(j/7r)r(i/)e*,2r 

(h+ioy 


(4.4.3) 


Hence,  for  the 

Sharp  right-angled  corner  (Figure  4.4a):  v  =  |,  and 


/  \  2  2iri 

w(z)  —  z*e  3  , 


1(h) 


-2r(fyt 

y/Z  (k  +  *0)s 


(4.4.4) 


The  Corcos  model  (3.2.4)  then  yields 


(U/8*)$(x,u)  2{2(q  A.  4 

mt  "  M‘  ( w) cos  (s(®  -  h  sm‘  *  s 


a’M./U) 


[al  +  (u6,/U)2]i 


3  5 


oo,  (4.4.5) 
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u8*IU 


Figure  4.4.  Predicted  noisespectrum((7/<5»)$(x,u;)/[ao(p0v2)2.M3(L<$*/|x|2)sin3  ^>cos2(|(0— 7r))] 
for  square-edged  and  rounded  right-angled  wedges;  R  is  the  mean  radius  of  curvature 
of  the  rounded  edge. 
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where  a'g  =  0.033.  This  approximation  is  strictly  applicable  when  both  faces  of  the  wedge 
are  large  compared  to  the  acoustic  wavelength. 

Rounded  right-angled  corner:  An  estimate  of  the  influence  of  rounding  can  be  obtained 
from  the  case  illustrated  in  Figure  4.4b,  for  which 

(4-4-6) 

where  R  is  the  mean  radius  of  curvature  of  the  edge,  and  C(z)  and  F(kiR)  are  defined  by 

1.18  —  0.126£  +  2.3£2  +  0.5£3  \  _  f 

(1  +  0.16f  +  O.lf2  +  £3)(1  +  0/  *  U  ’ 

(4.4.7) 

The  first  of  equations  (4.4.7)  maps  the  upper  half  of  the  C-plane  onto  the  profile  in  the 
plane  z-x i  +  ix2  of  the  fluid  region  bounded  by  the  rounded  wedge,  with  the  real  (-axis 
corresponding  to  the  wedge  profile.  The  acoustic  pressure  spectrum  $(x,w)  for  this  case  is 
obtained  by  multiplying  the  sharp-cornered  expression  (4.4.5)  by  aR/Uc). 

Figure  4.4c  compares  the  edge  noise  spectra  for  these  two  cases.  For  the  rounded  edge 
$(x,u;)  does  not  decrease  exponentially  fast  at  high  frequencies,  but  $(x,w)  ~  0( l/u1*). 
This  is  because  the  edge  rounding  in  the  model  of  Figure  4.4b  has  a  continuous  slope,  but 
the  curvature  is  discontinuous  at  each  end  of  the  rounded  sections. 

§4.5  Influence  of  mean  flow  direction 

When  the  mean  flow  over  the  trailing  edge  is  inclined  at  an  angle  j3  (<  90°)  to  the  aq-axis 
(see  Figure  4.2)  the  edge  noise  spectrum  continues  to  be  given  by  the  general  formula  (4.2.5). 
To  derive  the  analog  of  equation  (4.2.7)  for  the  acoustic  pressure  spectrum  we  must  make 
use  of  the  following  corrected  form  of  the  relation  (3.9.4) 

(Ps(ki,k3,u)p*(k[,  k3,u'))  «  — -5(u-u')6{ki  -  k[)P'{kuk3,u),  L^>5.  (4.5.1) 

Z7T  COS  fj 

where  L  is  the  transverse  width  of  the  turbulent  stream  (so  that  the  wetted  edge  of  the  airfoil 
has  length  L/  cos/3),  and  P'(ki,  k3,u)  here  denotes  the  blocked  pressure  spectrum  defined 
with  respect  to  the  x-coordinate  system.  This  system  is  inclined  at  angle  (3  to  the  mean 
flow.  Let  P(Ki,  K3,u>)  denote  the  wall  pressure  spectrum  defined  relative  to  a  coordinate 
system  in  which  the  K\- axis  is  parallel  to  the  mean  flow,  then 

P^ij^ocosr/qw)  «  P(fciCOs/?,  Aqsin/?,(j),  M  — >  0,  (4.5.2) 


z  _  —8 i 
R  37r 


(c»  +  a  +  0*).  r(hR) 
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and  (4.2.7)  becomes 


$(x,w) 


cuL  sin2(fl/2)  sin^  r°° 
27T2C0COS/3|x|2  J-  oo 


/OO 

\X(ki)\2P(ki  cos /?,  ki  sin /?,  u)  dki,  |x| 

-OO 


uL  sin2  (6/2)  sin  ip  (  v  \  r°° 

2-k2c0  cos2  /3|x|2  \Uc  cos  ft)  J- oo 


(4.5.3) 


The  Corcos  approximation  (4.3.7)  then  reduces  to 


(£7/<S,)$(x,  w) 

(PoV2)2 


(L5A  sin  sin2  (|)  /  u  \  (  u 

llxl2)  cos/?  \TrUccosp)  \Uccos(3 


_ aJj£8*/U) _ 

[oi2  +  M./C7)2]l[  1  +  (1-4)2  tan2  /?] 


,  a0  =  0.035.  (4.5.4) 
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5.  BLADES  AND  BLADE  TIPS 

The  efficiency  of  aerodynamic  sound  generation  by  lifting  surfaces  of  compact  chord  is  the 
same  as  for  an  acoustically  rigid,  compact  body  (a  M3).  The  amplitude  of  the  self-noise 
generated  by  a  short  section  of  the  trailing  edge  generally  depends  on  the  orientation  of  the 
section  relative  to  the  mean  flow  (§4.5).  The  importance  of  this  can  be  worked  out  exactly 
for  surfaces  of  elliptic  planform,  and  approximately  for  any  slowly  varying  planform.  In 
particular,  the  influence  of  blade  tip  contributions  can  be  estimated. 

§5.1  Green’s  functions 

We  consider  Green’s  functions  for  a  sequence  of  three  thin  airfoils  of  increasingly  complex 
planform,  but  of  acoustically  compact  chord. 

The  strip  airfoil: 

Let  the  airfoil  have  constant  chord  2a  and  infinite  span  (Figure  5.1a).  Take  coordinate  axes 
as  in  the  Figure,  with  the  origin  on  the  midspan  and  the  X\,  x2,  £3 -axes  respectively  in  the 
direction  of  the  mean  flow,  the  mean  lift  and  the  span.  The  time  and  frequency  domain 
compact  Green’s  functions  are  defined  as  in  (2.4.13)  and  (2.4.14)  with 

Xi  =  xi,  X2-  Re{  -  i\f£  -  a2],  X3  =  £3,  where  z  =  xi+  ix2  (5.1.1) 

Airfoil  with  slowly  varying  chord: 

The  definition  (5.1.1)  is  applicable  also  when  the  chord  varies  slowly  with  £3,  i.e.,  when 

y  .  da 

a  =  a(x3),  - —  <  1. 

dx  3 

Elliptic  airfoil: 

Let  the  planform  be  the  ellipse  (Figure  5.1b) 

/(£  1 ,  £3)  =  -2+-f  —  1  =  0,  (5.1.2) 

CL  C 

then  X\ 

¥>*(*) 

^(£) 
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£1,  X2—x2-  <p*(x),  X3  =  £3, 

'.i- _  * 
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(1  +  a2  +  y?)dix 


-x2 
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Figure  5.1.  (a)  Strip  airfoil  of  chord  2a;  (b)  Elliptic  airfoil  of  aspect  ratio  c/a. 
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and  A(x)  is  the  positive  root  of  [45] 


xj  ,4,  4 

a2  +  A  A  c2  +  A 


=  1. 


(5.1.4) 


Near  the  airfoil  A(x)  — >  +0,  and 

f{x  i,£3)H(/(:ri,:E3)) 


x\/aA  +  x§/c4 


+ 


Xn 


Xo 


X2 

*(o)  [VK 


i  +  ^Ci  +  A') 

,/A  2a  V  qV 


where  H  is  the  Heaviside  step  function. 


(5.1.5) 

(5.1.6) 


§5.2  Trailing  edge  self- noise 

Consider  low  Mach  number  mean  flow  at  speed  U  in  the  xi-direction.  The  self-noise 
generated  at  the  trailing  edge  can  be  expressed  in  the  form  (4.2.4).  In  the  acoustic  far  field 
the  analog  of  (4.2.5)  becomes 

p(x,  t )  *  f [°°  0Jj(k)l(k)ps(k,  a,)e-MH*l/* )dkldkzdu,  |x|  oo.  (5.2.1) 

o7TC0|x[  J  J — oo 

In  this  formula  k  =  (ku0,k3),  7  is  defined  as  in  (4.1.4),  0  =  arccos(x2/|x|)  is  the  angle 
between  the  airfoil  normal  (x2-axis)  and  the  radiation  direction,  and 


J(k)  =  j  [x2(y)]y^yi+{k3-KoCOS^y3}dS(y),  (5.2.2) 

where  $  is  the  angle  between  the  spanwise  direction  (x3-axis)  and  x,  and  the  integration  is 
over  the  upper  surface  (x2  =  +0)  of  the  airfoil.  The  function  AT2(y)  exhibits  the  following 
discontinuous  behavior  across  the  airfoils: 


[A^2(y)]_  =  X2{yi, -4-0, 2/3)  -  X2{yu  -0, y3) 


(5.2.3) 


=  2^a2  —  x2, 

=  2\Ja{xz)2  —  x2, 

2a  L  x2!  x| 

tf(0)V  a2  c2’ 


strip  airfoil; 
slowly  varying  chord; 

elliptic  airfoil. 


(5.2.4) 
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When  the  slowly  varying  chord  approximation  is  applied  to  the  elliptic  airfoil  we  make  the 
substitution 

a(x3)  ->  c 

in  which  case 

[*2(y)]*  ~  2 a 

which  agrees  with  the  exact  expression  in  (5.2.4)  provided  (0)  ~  1.  Figure  5.2  shows  that 
this  is  an  excellent  approximation  when  the  aspect  ratio  a  =  c/a  of  the  elliptic  airfoil  exceeds 
about  5. 

To  evaluate  the  integral  (5.2.2)  only  the  contribution  from  the  vicinity  of  the  trailing  edge 
xi  =  a(x 3)  must  be  retained.  For  the  strip  airfoil,  and  turbulence  length  scales  that  are 
much  smaller  than  the  airfoil  chord  ( ka  1),  we  write 

X(k)  =  2  j  \Ja?  -  y\  e<{fc>w+(fc»-'t»CM^w}dS(y) 

ro. 

~  4irV2a6(k3  -  K0cosip)  /  ^ a  -  y[  eiklVldyi  =  -27r\Z27ra<5(fc3  -  k0  cos  ip)  — 

J —oo  rC^ 

(5.2.5) 

Using  this  in  (5.2.1)  and  calculating  the  acoustic  pressure  frequency  spectrum  $(x,a>)  for 
the  wetted  section  of  the  trailing  edge  of  span  L  shown  in  Figure  5.1a,  we  find  (using  (3.9.4)), 
for  the 


Strip  airfoil: 


$(x,w)  ~ 


aL  oj2  cos2  0  r°°  P(fci,K0 cos  dk\ 
8|x|2  c2  J-o o  |fci| 


n  aLM2  cos2  Q^pnioj)  U  _ ,  ,  ,com 

0.027 - j— tz — M  =  —  <  1,  ua/U  »  1,  |x|  ->•  oo,  (5.2.6) 


where  the  Corcos  formula  (3.2.3)  has  been  used  to  obtain  the  second  line. 


The  integral  (5.2.2)  defining  J(k)  for  the  elliptic  airfoil  is  evaluated  by  first  expanding 
[^2(y)]l  about  the  trailing  edge  y\  =  aJl  —  yf/c2,  as  in  (5.2.5),  by  writing 


Section  5.2 


69 


Blades  and  Blade  Tips 


Report  No.  AM  98-001 


Boston  University,  College  of  Engineering 


0  10  20  30  40  50 

aspect  ratio  c/a 


Figure  5.2.  Dependence  of  'I'(O)  and  1/1F2(0)  on  aspect  ratio  c/a. 
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and  performing  the  ^-integration  over  the  range  -oo  <  3/1  <  a\Jl  -  y\/c2.  Let  us  suppose, 
as  indicated  in  Figure  5.1b,  that  a  section  of  the  trailing  edge  of  length  L  inboard  of  the 
blade  tip  is  wetted  by  the  flow,  and  examine  the  contribution  to  the  radiation  from  the  tip 
region.  Then 

z(k>  *  ?tfP  £1 0  ■  exp  {*  (haf~^+ ~  k»c°s^>!'3)  }  *»■ 

When  the  turbulence  scales  are  much  smaller  than  the  airfoil  chord  {ka  »  1),  the  principal 
contributions  to  this  integral  are  from  the  neighborhood  of  that  value  of  y3  where  the  phase 
of  the  exponential  is  stationary.  Geometrically,  this  corresponds  to  the  condition  that  the 
vector  (Jfci,  0,  A;3  -  k0  cos  ip)  be  parallel  to  the  normal  to  the  elliptic  trailing  edge,  and  yields 

-27to2c  exp  +  (fc3  -  k0  cos  ip)2 c2}  f  k3  -  k0  cos  ip  q(l  —  L/c)  \ 

~  ^(0)[A:2a2  +  (k3  —  K0cosip)2c2]  \  h  \[Lc^2  -  L/c / 

(5.2.7) 

H  is  the  Heaviside  step  function,  whose  argument  is  positive  provided  the  point  on  the  trailing 
edge  where  (fci,  0,  k3  -  k0 cos  ip)  is  normal  to  the  edge  lies  within  the  interval  c-  L  <  x3  <  c 
wetted  by  the  flow.  It  may  be  verified  that  this  result  reduces  to  (5.2.5)  when  c/a  00. 

Hence,  evaluating  the  acoustic  pressure  frequency  spectrum  in  the  usual  way  (making  use 
of  (3.1.3)),  and  noting  that  k0  <C  k,  we  find  for  the 


Elliptic  airfoil: 

a4c2w2  cos2  0  f°°  fc2P(k,  w) 


$(x,w) 


8$2(0)  C; 


os  fc)  f°°  k 

WJ-ooWl 


H  |  fa  _  o(l  -  L/c) 


{k\o?  +  kjc2)2  \kx  JTcfi^L/c, 


dkidkz,  |x|  ->•  00. 

(5.2.8) 


Using  the  Corcos  approximation  (3.2.3),  this  result  can  be  expressed  in  the  form 

aLM 2  cos2  Q$pp(a;) 


K 


4(x,u>)  «  *„(x,w)k(|,|)  ,  *„(*.«>)  s  0.027 
/  c  c  \  _  1  c  c  r00  (1  +  /x2)d/x 

Vo’lJ  “  mouiJv 


tf2(0)  a  L  J„0  (1  +  (1.4)2/x2)(l  +  cV/a2)2  1 


a  (c/L  —  1) 

Vo  =  -- 
c 


y/2c/L-l 


(5.2.9) 


Thus,  JC(c/a,  c/L)  represents  the  ratio  $(x,  w)/$0(x,w)  of  the  edge  noise  spectrum  to  that 
of  the  sound  generated  by  a  trailing  edge  section  of  span  L  of  a  strip  airfoil  of  chord  2 a.  This 
is  plotted  in  Figure  5.3  as  a  function  of  the  fractional  span  of  the  wetted  tip  L/ c  for  three 
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different  aspect  ratios  c/a.  The  contribution  from  the  tip  region  is  6  dB  or  more  below  that 
from  an  inboard  section  of  the  airfoil  of  chord  2 a  and  equal  span  when  L/c  <0.1.  According 
to  (5.2.9),  K.  ~  tt/4#2(0)  »  7r/4  provided  a  <  c  and  L  =  c„  i.e.,  the  maximum  edge  noise 
from  the  elliptic  airfoil  is  always  about  1  dB  smaller  than  that  from  an  airfoil  of  uniform 
chord  a  of  the  same  span.  Note  that,  since  tf2(0)  «  1  for  c/a  >  5,  these  conclusions  remain 
true  when  the  calculation  is  based  on  the  slowly  varying  chord  approximation. 
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Figure  5.3.  Edge  noise  produced  by  a  tip  region  of  length  L  of  an  elliptic  airfoil  of  maximum 
chord  2 a  relative  to  that  produced  by  the  same  flow  over  a  spanwise  section  of 
length  L  of  a  strip  airfoil  of  chord  2o. 
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6.  SLOTS  AND  FLAP  SIDE-EDGES 

The  side  edges  of  part  span  trailing  edge  flaps  are  important  contributors  to  trailing  edge 
noise  [46  -  51].  In  typical  underwater  applications  the  streamwise  “slot”  between  an  un¬ 
deployed  flap  and  the  neighboring  fixed  section  of  a  trailing  edge  is  large  enough  to  be  an 
important  additional  source  of  noise.  This  case  can  be  handled  analytically  by  modelling 
the  trailing  edge  region  by  a  rigid  half-plane  (as  for  conventional  trailing  edge  noise,  §4)  and 
the  slot  as  a  rectangular  cut-out,  as  indicated  in  Figure  6.1a.  In  a  first  approximation  the 
influence  of  flap  deployment  on  the  generation  of  sound  can  be  estimated  by  assuming  the 
width  2s  of  the  slot  to  increase  with  distance  towards  the  trailing  edge. 

§6.1  Green’s  function:  non-compact  airfoil  chord 

The  slot  is  assumed  to  be  acoustically  compact  and  to  occupy  the  region  -l  <  x\  <  0,  x2  = 
0,  |z3|  <  s  of  the  rigid  half-plane  xx  <  0,  x2  =  0.  For  an  observer  at  x  in  the  acoustic  far 
field,  the  compact  Green’s  function  is  written 

G(x,  y;  u)  =  G0(x,  y ;  u)  +  <?i(x,  y;  u)  +  Gs(x,  y;  w),  (6.1.1) 

where  G0,  G i  are  defined  as  in  (4.1.2)  for  the  rigid  half-plane,  and  Gs  is  the  correction 
produced  by  the  slot. 

The  functional  form  of  Gs  depends  on  the  minimum  distance  A  of  the  source  point  y  from 
the  slot,  i.e.,  on  the  magnitude  of  A/s: 


For  A  ~  s: 


Gs(x,y;w) 


For  A  »  s: 


sin  2  V’  sin  (0/2)  y/ K0telKo  M 
7T\/27ri|x| 


/%i)Re  <  In 


{■ 


,  Z  =  2/3  +  iy  2; 

(6.1.2) 


<?s(x,y;o;) 


2sin2  ^sin(0/2)\//c^£el'Colxl  f°  /?(£) 


7r2\/27ri|x| 


_  — sgn(y2)  sinz  if)  sin(0/2)\/K^el'c<’lxl 

87r3-\/27ri|x| 

f0  ^  [CO  ei[ki+km+km+'tM](ikdkidk3  ftt  1 

+  h  +  «0) ;  {"} 
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slot 


Figure  6.1.  (a)  Model  of  a  trailing  edge,  streamwise  slot; 

(b)  View  of  the  slotted  half-plane  from  downstream. 
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where 


(yi,y2)  =  r(cos0o,sin0o), 


a  = 


2r|fl  -  sin  - 


\  R{R  +  Ri)  2  ’ 

#  =  \/(y i  —  f)2  +  2/2  +  yi)  =  \Ar  ■**  l£l)2  2/3) 


and  7  and  R0  are  defined  as  in  (4.1.4).  The  derivation  of  these  formulae  is  discussed  in  §6.3. 

The  angles  ip  and  6  specify  the  far  field  observer  direction  x,  as  in  Figure  4.2.  The 
dimensionless  function  /3(xi)  defines  a  rectilinear  source  distribution  to  which  the  slot  is 
acoustically  equivalent  when  irradiated  by  a  (reciprocal)  point  source  placed  at  x.  This 
function  is  tabulated  in  Table  6.1  and  plotted  in  Figure  6.2  for  four  different  values  of  the 
slot  aspect  ratio  2 s/£. 


2s/t  — » 

0.025 

0.05 

0.075 

0.1 

Xi/l 

P(xi) 

-1.0 

0.310 

0.360 

0.390 

0.420 

-0.9 

0.234 

0.282 

0.310 

0.344 

-0.8 

0.204 

0.240 

0.268 

0.291 

-0.7 

0.184 

0.216 

0.239 

0.259 

-0.6 

0.167 

0.195 

0.216 

0.234 

-0.5 

0.151 

0.177 

0.196 

0.212 

-0.4 

0.136 

0.159 

0.176 

0.190 

-0.3 

0.119 

0.139 

0.155 

0.168 

-0.2 

0.100 

0.118 

0.132 

0.143 

-0.1 

0.077 

0.090 

0.103 

0.113 

0.0 

0.000 

0.000 

0.000 

0.000 

Table  6.1 


The  approximation  (6.1.2)  is  applicable  for  source  distributions  in  the  immediate  neigh¬ 
borhood  of  the  slot  whose  length  scale  in  the  ^-direction  is  comparable  to  the  slot  width 
2s,  whereas  (6.1.3)  should  be  used  when  the  bulk  of  the  source  distribution  is  spread  over 
a  much  larger  spanwise  domain  which,  nevertheless,  is  well  within  the  acoustic  near  field  of 
the  slot. 
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X\lt 


Figure  6.2.  Dependence  of  the  Green’s  function  source  strength  (3(x i) 
on  the  slot  aspect  ratio  2 s/t. 


Section  6.1 


77 


Slots  and  Flap  Side-Edges 


Report  No.  AM  98-001 


Boston  University,  College  of  Engineering 


§6.2  Green’s  function:  compact  chord 

For  a  thin  airfoil  whose  chord  at  the  slot  is  2a  »  l  (Figure  6.3)  the  compact  Green’s 
function  in  the  absence  of  the  slot  is  given  by  (2.4.14),  when  the  coordinate  origin  is  at 
midchord.  When  the  origin  is  shifted  to  the  trailing  edge,  with  the  a; i -axis  along  the  axis  of 
symmetry  of  the  slot,  the  analogs  of  equations  (6.1.2),  (6.1.3)  are  found  to  be: 


For  A  ~  s: 


i  cos  0\/2/c^a£etK°lxl 
Gs(x,y;w)« —  v - 


47t|x| 


L 

2  +  ^-1 

l 

s  Vs2 

,  z  =  y3  +  m ;  (6.2.1) 


For  A  »  s: 

Gs(x,y;w) 


— tcos 


bi«o|x| 


27t2|x| 


/: 


i  R(Z,  y) 


— ) 
+  a2/ 


zsgn(y2)  cos  0^2/c2a(?e,'Colxl 
~  327t3|x| 

,o  roo  em+km+km)-M n  0, 

TJmmmFThT®)  ;  (62'2) 

where  0  =  cos-1(x2/|x|)  is  the  angle  between  the  radiation  direction  and  the  upward  normal 
to  the  airfoil  (the  x2-axis). 


§6.3  Derivation  of  Green’s  function 

We  use  the  method  of  [49]  making  use  of  the  reciprocal  theorem,  in  which  G(x,  y;  u>)  = 
G(y,x]u)  is  first  calculated  as  a  function  of  y  in  the  neighborhood  of  the  slot  when  a  unit 
point  source  is  placed  at  the  far  field  point  x  (defined  as  in  (2.4.2)).  Details  are  given  below 
for  an  airfoil  of  non-compact  chord,  but  the  procedure  is  readily  modified  for  the  case  of 
compact  chord  discussed  in  §6.2. 

The  length  i  of  the  slot  is  assumed  to  be  small  compared  to  the  acoustic  wavelength.  Thus,  when  scattering 
by  the  slot  is  ignored,  the  potential  in  the  neighborhood  of  the  slot  produced  by  an  acoustic  source  at  the 
far  field  point  x  is  given  by  (4.1.1).  In  the  immediate  vicinity  of  the  slot  the  motion  induced  through  the 
slot  is  locally  two-dimensional  and  can  be  approximated  by 

G  =  A  +  B Rej  ln[z/s  +  \/z2/52  -  l]j,  z  =  r/3  +  iy2,  (6.3.1) 

where  A  =  ,4(x;w)  is  independent  of  y,  and  B  =  B(x,yi;w)  determines  the  local  volume  flux  through  the 
slot  (which  varies  with  y\  along  its  length)  induced  by  the  incident  potential  (4.1.1). 
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Section  6.3 


Figure  6.3.  Trailing  edge  slot  in  airfoil  of  compact  chord. 
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At  large  distances  (»  s)  from  the  slot  the  influence  of  the  induced  motion  through  the  slot  is  equivalent  to 
that  produced  by  a  line  source  of  variable  strength  <7(11),  say,  per  unit  length  distributed  over  -i  <  xx  <  0 
on  the  upper  surface  x2  =  +0,  together  with  that  produced  by  an  equal  and  opposite  distribution  (-a(xi) 
per  unit  length)  on  the  underside  x2  =  -0.  At  distances  from  the  slot  that  greatly  exceed  s  but  are  well 
within  an  acoustic  wavelength,  the  net  field  can  therefore  be  cast  in  the  Macdonald’s  form  [52,  53  (§8.4.2)] 


G  =  Go  +  Gi  — 


tan  1 


( 


y/l  +  Ot2  J 


<*e, 


(6.3.2) 


where  R  and  a  are  defined  as  in  (6.1.3). 


It  follows  from  a  comparison  of  (6.3.1)  and  (6.3.2)  that,  since  K0yz  4C  1, 


A  =  G0 


47t|x| 


(6.3.3) 


To  determine  B  and  cr  we  introduce  a  length  h  which  is  assumed  to  be  large  compared  to  s,  yet  small 
compared  to  the  length  scale  of  variation  of  B  and  a  along  the  axis  of  the  slot.  Then,  for  each  fixed  value  of 
j/x  in  the  interval  -i  <  yi  <  0  we  find  from  (6.3.1),  (6.3.2) 


G 


A  +  B  In 


\z\/s  -»  00 


Go  + <3,  +  £<!!!>  1„ 


7T 


-  +  —  I  — 

h  2n  J_(  \yi 


£1 


4 

1 - tan 

7 r 


M=)l 

\\/l  +  a2/. 


^  2 nj-t  iw-er 


\z\/h 


a=vi^ 


2£yi 


^Klyi  - -el  -  yi  - -e)  ’ 


(6.3.4) 


where  the  principal  value  notation  implies  that  the  interval  j/i  —  h/2  <  f  <  yi  +  h/2  is  to  be  excluded  in  the 
final  integral.  Hence  B  =  cr/ir,  and  by  setting 


-sin^sin(fl/2)%/^fei*°lxl 

*V2ri\x\ 


m, 


(6.3.5) 


it  follows  that  P  is  the  solution  of  the  integral  equation 


_  1  f°  P(rp 

K-*?l  2  -  77| 


-t  <  Z  <  0.  (6.3.6) 


This  can  be  solved  by  collocation  (Figure  6.2);  the  solution  is  stable  provided  h/2s  >  1. 


§6.4  Boundary  layer  generated  slot  noise 

The  sound  generated  by  turbulent  flow  in  the  Xi-direction  over  the  slot  can  be  calculated 
by  the  method  of  §4.2,  by  using  the  formulae  (4.2.3)  and  (4.2.4)  to  express  the  radiated  sound 
in  terms  of  the  blocked  pressure  ps.  Because  the  dominant  boundary  layer  pressures  have 
surface  wavenumbers  (ki,  k3)  in  the  convective  region  centered  on  ( u)/Uc ,  0),  the  characteristic 
length  scale  of  the  source  field  in  the  spanwise  direction  greatly  exceeds  the  slot  width  2s.  It 
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is  therefore  appropriate  to  calculate  the  slot  noise  by  using  one  of  the  “outer”  forms  (6.1.3), 


(6.2.2)  of  Gs. 

For  an  airfoil  of  non-compact  chord,  equations  (6.1.3)  and  (4.2.4)  yield  the  far  field  acoustic 
pressure  in  the  form 


K  sini^inff/2)  rr°°  dikdu  r°  Mmel{  (,-<f ^^kl) 

’  s/2i r?x  Jj- OO  j-e  v  ' 

1  (6.4.1) 


where  [t]-t-  |x|/c0  is  the  retarded  time.  Forming  the  acoustic  pressure  frequency  spectrum 
in  the  usual  way  (using  (3.1.3))  we  find 


*(«,«)  *  ut*u*'£m  CHk,w) 


drk 


ui  sin  tp  sin2  (9/2) 


7TC0|x|2 


'*PpM 


w  >0, 

(6.4.2) 


where  $pp(u;)  is  the  blocked  pressure  frequency  spectrum. 

Using  the  Chase  [27]  approximation  (3.2.4)  this  result  can  be  written 


(£//5»)$(x,  u)  i2  ...  2 

[»..;)■  *“g  *"»■• 


(5) 


■7»M*/u)2 


where 


(6.4.3) 


(6.4.4) 


is  a  dimensionless  function  of  the  frequency.  The  asymptotic  formula  for  the  behavior  of  the 
error  function  as  uj£/Uc  -¥  oo  [54]  supplies 

r°  mdd 


HH 


~  Foo  = 


ttUc 
20  Ul 


i: 


L Ot 

~uc 


oo. 


(6.4.5) 


Thus,  $(x,  u)  decreases  like  ui  1  at  high  frequencies,  whereas  the  trailing  edge  noise  generated 
by  the  same  flow  decreases  at  least  as  fast  as  u;-2  (see  (4.3.7)). 

The  situation  is  illustrated  in  Figure  6.4a  for  a  slot  of  aspect  ratio  2 sfl  —  0.05,  and  for 
i/8  =  10,  where  8  ~  85*  is  the  boundary  layer  thickness.  The  dashed  curve  is  the  spectrum 
(C//5*)$(x,  w)/[a0(p0u2)2M(!$*/|x|2)  sin xjj  sin2 (6/2)]  (dB) 
of  the  sound  generated  by  a  section  of  the  trailing  edge  whose  spanwise  length  equals  the  slot 
length!  The  “slot  noise”  is  the  spectrum  {U/8^(x,u)/[a0(p0vl)2M(£8t/\x\2)  sin  ip  sin2  {d/2)] 
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Figure  6.4.  (a)  Comparison  of  slot  noise  and  the  trailing  edge  noise  generated  by  a  section 
of  the  edge  of  span  l.  (b)  Dependence  of  slot  noise  on  i/8. 
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determined  by  (6.4.3)  for  Uc  -  0.7 U.  At  very  high  frequencies  the  amplitude  of  the  slot  noise 
fluctuations  (caused  by  the  frequency  dependence  of  F(u))  decrease,  and  the  high  frequency 
limit  shown  in  the  figure  is  the  prediction  of  (6.4.3)  when  ^(w)  is  replaced  by 

Evidently  the  high  frequency  approximation  provides  a  good  overall  representation  of  the 
strength  of  the  slot  noise  for  all  frequencies,  and  it  has  been  used  in  Figure  6.4b  to  illustrate 
the  dependence  of  slot  noise  on  i/5. 

§6.5  Side-edge  flap  noise 

To  investigate  the  sound  generated  by  the  side-edge  vortex  formed  when  a  part-span  flap  is 
deployed  at  a  finite  angle  of  attack  (Figure  6.5a),  it  is  necessary  to  modify  the  slot  geometry. 
The  distance  between  the  adjacent  sides  of  the  flap  and  the  undeflected  trailing  edge  now 
increases  as  the  trailing  edges  are  approached;  in  a  first  approximation,  the  influence  of 
flap  deflection  at  the  side-edge  can  therefore  be  modeled  by  taking  the  slot  width  2 s  for  an 
undeflected  flap  to  vary  along  the  slot.  The  simplest  case  is  that  in  which  flap  deployment 
consists  of  a  simple  rotation  through  a  small  angle  d  about  the  flap  leading  edge,  as  illustrated 
in  Figure  6.5a.  In  this  case  flap  deployment  is  approximated  by  the  method  of  §§6.1  -  6.3 
by  assuming  the  slot  semi-width  s  =  s(y i)  to  increase  with  yi  (Figure  6.5b)  according  to 

s(yi)  =  s0  [l  +  {y\jt  +  1)  (0.  +  (^£/2s0)2  -  l)] ,  —i  <  y\  <  0,  (6.5.1) 

where  2s0  is  the  slot  width  in  the  undeflected  state  (i9  =  0). 

The  compact  Green’s  function  is  given  by  the  formulae  of  §§6.1,  6.2,  except  that  explicit 
account  is  now  be  taken  of  the  dependence  s  on  j/j.  Equation  (6.3.6)  for  the  source  strength  ft 
is  unchanged  provided  s  is  replaced  by  s(£).  The  effect  of  this  modification  on  the  functional 
form  of  P(xi)  is  illustrated  in  Figure  6.6  when  2 s0/t  —  0.05  for  $  =  0°,  5°,  10°.  The 
magnitude  of  (d{x\ )  is  seen  to  be  determined  principally  by  the  local  slot  width  2s:  fl{x{) 
differs  significantly  from  its  value  in  the  undeflected  state  only  towards  the  trailing  edge. 

To  estimate  the  side-edge  noise  generated  by  the  flap  lift  vortex  T ,  which  begins  to  form 
near  the  forward  part  of  the  side-edge  (Figure  6.7),  assume  the  mean  position  of  the  vortex 
lies  along  the  line  X\  >  —i,  x^  «  $(yi  +  £),  x$  =  s0,  where  1  and  the  xi-axis  is  taken 
to  be  parallel  to  the  flap  edge.  For  small  x)  the  component  of  the  unsteady  flow  velocity 
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Figure  6.5.  Modelling  flap  deployment  by  a  side-edge  slot  of  variable  width. 
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Figure  6.6.  Dependence  of  f3(x\ )  on  when  2s0/i  =  0.05. 
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normal  to  the  vortex  axis  is  taken  in  the  form 

vn  =  vn(t  -  Xl/Uc)  =  (0,  v2,  vs)  =  f°°  vn(u)eMxi/U'~t)fo  (6-5-2) 

where  Uc  is  a  suitable  convection  velocity  associated  with  the  propagation  of  disturbances 
along  the  vortex  core.  Using  the  general  solution  (2.6.4)  (with  v  =  0  and  VjUj  =  0  on  S), 
with  Gs  given  by  (6.1.2)  for  an  airfoil  of  non-compact  chord,  the  sound  generated  by  the 
interaction  of  the  vortex  with  the  side-edge  is  found  to  be  given  by 

«  1, 

(6.5.3) 

where  [t]  =  t-  |x|/c0  is  the  retarded  time. 

Let  v„(t  —  Xi/Uc)  be  a  stationary  random  function,  and  let  $23 (w)  >  0  be  the  frequency 
spectrum  defined  by 

((v3(u>)  -  v2(w))  (v^w')  -  v2V)))  =  U25(tv  -  w')$23(u),  (6.5.4) 


_  p0rsin^sin(g/2)  p  r°  g(yQvggfaM  ~  V2{w)}eMxi/Uc  [t])dyi 
27r\/27ri|x|\/^  J-00  J-t  \/s(yi)(j/i  +1) 


where  the  asterisk  denotes  complex  conjugate.  The  side-edge  noise  frequency  spectrum  can 
then  be  written 

*(x,  w)  «  i23(w)s(u>eiUc)t  |x|  _>  oo,  u  >  0,  (6.5.5) 

where  Mc  =  Uc/c0,  and  the  dimensionless  efficiency  factor  £ (cv£/Uc)  is  given  by 


(6.5.6) 


£  =  ue  r° _ _ (65.6) 

Uc  J-t  [1  +  (yjl  +  1)(^1  +  (<M/2».)S  -  1)]  JUvi  +  t) 

£(u£/Uc)  turns  out  be  effectively  independent  of  the  flap  angle  d  (at  least  for  d  <  10°),  and 
is  plotted  in  Figure  6.7  for  2 s0/i  =  0.05. 

In  a  rough  approximation  the  vortex  strength  r  ~  Uid,  and 

.  McplU4i2  sin  ip  sin2(0/2)  / \  ,  . c,  .  .  (rk7\ 

$(x,u;)«— ^ (  57  1  $23{w)S(ue/Uc),  M-+00,  (6.5.7) 


wherein  id/2s0  is  the  ratio  of  flap  deflection  at  the  trailing  edge  to  the  undeflected  slot  width. 
Further  predictions  from  this  formula  require  a  quantitative  knowledge  of  the  frequency 
spectrum  $23 (w). 
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In  addition,  however,  the  side-edge  vortex  also  produces  sound  by  interacting  with  the 
trailing  edges  of  the  flap  and  the  neighboring  undeflected  airfoil,  determined  by  the  com¬ 
ponent  G\  of  the  composite  Green’s  function  (6.1.1)  (see  (4.1.2)).  This  can  be  significant 
only  at  lower  frequencies,  when  the  length  scale  of  the  associated  unsteady  hydrodynamic 
pressures  are  large  enough  to  affect  these  edges.  The  principal  source  term  in  (2.6.4)  may 
now  be  approximated  by  u>  A  v  =  (0,  -u3I\  0)%2)<%3),  produced  by  spanwise  oscillations 
of  the  vortex  core.  In  this  case  we  find 


p(x,  t ) 


—p0Ty/M~csm*  xpsm(Q/2) 

2\/27r|x| 


u3(t  -  |x|/c0), 


00, 


(6.5.8) 


where  u3  is  defined  in  (6.5.2). 

Hence,  forming  the  acoustic  pressure  frequency  spectrum  as  before: 

Mcp20U2T2  sin  ip  sin2 (9/2) 


$(x,w) 


47r2|x|: 


$3(u>) 


$2Mcp2f74^2  sin  xp  sin2  {6/2) 
47t2|x|2 


$3(u;), 


oo, 


(6.5.9) 


where  (u3(u)ul(u'))  =  U25(u>  -  w')<h3(a;).  For  small  deflections  d  it  appears  that  this  will 
be  small  compared  to  the  side-edge  generated  sound  (6.5.7). 

The  results  (6.5.7)  and  (6.5.9)  represent  the  sound  produced  by  perturbed  motions  of 
the  vortex  T  without  taking  account  of  contributions  from  any  other  vorticity  that  may  be 
produced  at  the  surfaces  of  the  flap  and  airfoil  by  that  motion.  The  importance  of  these 
additional  sources  of  sound  is  not  known,  and  can  only  be  determined  from  a  more  detailed 
knowledge  of  the  hydrodynamic  flow  in  the  neighborhood  of  the  side  edge. 
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PRINCIPAL  SYMBOLS 


S  boundary  layer  thickness 

St  displacement  thickness 

r]  shear  coefficient  of  viscosity 

k0  acoustic  wavenumber  u/c0 

$  incompressible  velocity  potential  satisfying  (2.4.7) 

n  \/MM 

v  kinematic  viscosity 

7Tjj  momentum  flux  tensor 

$(x,  u)  acoustic  pressure  frequency  spectrum 

$pp(uj)  wall  pressure  frequency  spectrum 

p  fluid  density 

p0  mean  fluid  density 

crij  viscous  stress  tensor 

uj  vorticity  curl  v 

u!  radian  frequency 

B  total  enthalpy  w  +^v7 

c  speed  of  sound 

c0  mean  sound  speed 

D/Dt  d/dt  +  Vjdf  dxj 

Di  drag  force 

Fi  surface  force  applied  to  fluid 

k  planar  wavenumber  (ki,  0,^3) 

k  m 

K0  vacuum  bending  wavenumber 

M  Mach  number 

p  pressure 

p0  mean  pressure 

ps  wall  pressure 

Pij  compressive  stress  tensor 

P(k,u;)  wall  pressure  spectrum 
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RH 

n(xi,x3,t) 

M) 

s 

t 

T 

Tij 

U 

Uc 

v 

V» 

w 

X 


y 

Yt 


elastic  wall  refection  coefficient 
wall  pressure  correlation  function 
polar  coordinates  in  the  rcia^-plane 
specific  entropy 
time 

temperature 
Lighthill  stress  tensor 
mean  velocity 
convection  velocity 
velocity 

friction  velocity 
specific  enthalpy 
(xi,x2,x3) 

( yi,  2/2, 2/3 ) 

Xi  -  <p*i(x) 

Vi  -  <p,-(y) 
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